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ABSTRACT 


The drag on a body with a self-contained electromagnet are meas- 
ured with the body is placed in a low density argon plasma stream. 

Tests are conducted in a steady flow of plasma produced in a direct 
current arc and expanded to supersonic speeds. Measurements are re- 
stricted to one flow condition where the influence of viscosity and 
the Hall effect are both important. 

From experimental measurements of the total drag on the body, the 
Lorentz force on the magnet and the pressure force on the body, the 
viscous forces on the body may be deduced. Empirical relations for 
the data are then generated. It is found that the total drag and 
Lorentz force increase with increasing magnetic field but that the 
pressure and viscous forces decrease. A surprising result is that 
pressure and density increase at the stagnation point with magnetic 
field contrary to previous predictions. This indicates that when Hall 
effect is important, the assumption of constant density in the shock 
layer is no longer justified. The experimental results are compared to 
existing theoretical treatments of the problem. 

The flow field around the body is investigated spectroscopically, 
and the predicted increase in shock stand-off distance is observed 
photographically. An unexpected visible color change which occurs when 
the magnetic field is applied is also investigated. Spectroscopic meas- 
urements show that the electrons in the plasma are out of thermal equi- 
librium in the plasma and the flow is essentially non reacting. The in- 
fluence of these facts on other experimental measurements is considered. 
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CHAPTER I 
INTRODUCTION 


Entry of a hypersonic vehicle into the atmosphere is known to present 
severe problems caused by the shock-heated gas at the nose of the body. 

Not only is heat transfer to the body extreme but due to ionization in the 
gas, communications to the vehicle are difficult during part of the entry 
flight. 

Accompanying the interest in the problems of high speed entry is the 
proposal by Kantrowitz (Ref. 1) and by Resler and Sears (Ref. 2) to utilize 
the interaction of the field of an onboard magnet with the ionized flow in 
the forward region of the vehicle. Consideration of this interaction has 
shown that a variety of different effects will appear offering such possi- 
bilities as flight control based on the change in drag and lift on a 
vehicle with a magnetic field. Simultaneously, heat transfer to the body 
is reduced. 

Figure 1 depicts the changes in the hypersonic flow around a blunt 
body when the flow is ionized and a magnetic field is applied. Here the 
body is a simple hemisphere at zero angle of attack to the flow. The 
magnet is presumed dipole- like, with the axis aligned with the flow. Due 
to the flow of charged particles with velocity V across a magnetic field 
B, an electric field and consequently a current J is induced in the 
azimuthal direction in accordance with Ohm's Law. The interaction of this 
current with the magnetic field results in Lorentz force F on the fluid. 

It is seen that this force is given byF=JxB=aVxBxB where a is 
the conductivity of the gas and the force is directed to oppose the flow. 
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In a normal aerodynamic flow of this nature, the shock wave which forms 
in front of the blunt body is detached to permit passage of the flow around 
the body as dictated by continuity. Since the flow is slowed when the 
magnetic field is applied, the detachment distance must therefore increase. 
Likewise the viscous drag and heat transfer on the body must change since 
these quantities depend on the gradient of velocity normal to the body 
which is altered. A component of the force on the fluid lies along the 
direction of the freestream flow. This force partially supports the flow 
in the shock layers so that the pressure force on the body is relieved. 
Simultaneously the reaction to this component acts directly on the magnet, 
thus the drag force on the body is increased by this amount. 

In this explanation we have avoided the Hall effect by taking Ohm's 
Law in its simplest form, i.e., we have ignored the reduction in conduc- 
tivity in the direction of the electric field and the introduction of 
transverse currents, both due to the application of a magnetic field. A 
reduction in conductivity in the azimuthal direction will lead to a reduc- 
tion in any advantage gained by the application of a magnetic field. 
Furthermore, the additional Hall currents also interact with the magnetic 
field to produce a force in the azimuthal direction on the fluid and a 
torque on the magnet. 

We have only discussed a small part of the phenomena of magnetoaero- 
dynamics. An elementary consideration of magnetoaerodynamic lift could proceed 
in a similar manner with the body or magnet or. both taken at an angle of 
attack to the flow (see Ref. 3). In a practical flight situation a vehicle 
entering the earth's atmosphere passes through a regime of rarefied flow. 

Here electric forces on the body and their alteration by a magnetic field 
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are known to be important. For a discussion of these phenomena the reader 
is referred to Ref. 4. 

This study is a laboratory investigation of the magne togas dynamic 
flow around a blunt body at supersonic speeds. The fundamental aspects 
of magnetoaerodynamics are examined for a simple case. We restrict the 
study to an axisymmetric conducting body with either a hemispherical or 
flat end nose, at zero angle of attack. The axis of the solenoidal magnet 
is aligned with the flow. In addition to these restrictions, we concen- 
trate on a single continuum flow condition and examine in detail the mag- 
petoaerodynamic drag on the body and the alteration to the shock layer. 

Under these restrictions we have precluded an investigation of 
magnetoaerodynamic lift or electric forces on the body. Furthermore, in 
the laboratory it is expected that interaction of the flow and magnetic 
field leaves the magnetic field undistorted (Ref. 5). In other words, 
the appropriate parameter, the magnetic Reynolds number (to be discussed 
and evaluated later) is so small that distortion of the magnetic field 
lines is negligible. While we are primarily concerned with a laboratory 
problem, qualitative differences with a flight situation will be noted as 
they occur. 

The rest of this chapter contains a review of the literature per- 
taining to blunt body magnetoaerodynamic flows (summarized in Table I). 
Following this the equations of magnetogas dynamics are stated and a preliminary 
analysis of the problem is made. Chapter II presents a background for the 
laboratory investigation and Chapter III concerns the facility in which the 
experiments are to be performed. The experiments pertaining to magneto- 
aerodynamic drag are reported in Chapter IV. Similarity, Chapter V reports 
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the detailed examination of the flow field and measurement of magneto- 
aerodynamic shock standoff distance by optical methods. The results of 
Chapters IV and V are further discussed and comparison is made to theory 
in Chapter VI. Chapter VII contains a summary of the investigation and 
recommendations for future research. 


1.1 Literature Survey 

In 1958, following the initial suggestions of the possible uses of 
magne togas dynamics during entry, Bush (Ref. 6) treated the problem of 
magnetoaerodynamic hypersonic flow around a blunt body. The analysis is 
confined to the stagnation region and assumes an inviscid flow of constant 
density and constant conductivity behind the shock. Specifying the mag- 
netic field at the shock and substitution of suitable functions for the 
field and velocity permits reduction of the problem to a point where it 
may be handled numerically. Both the case of small and of vanishing 
magnetic Reynolds numbers are treated. Results are presented in terms of 
the magnetic interaction parameter defined using the magnetic field 

strength at the shock and tl^radius of the shock as the characteristic 

a l R s B s 

dimension. Thus, S = where a. is the conductivity behind the 

S 0 V X 

r 00 00 

shock and the subscript ©o denotes frees tream conditions. The results 
show an increase in shock standoff with field strength and a decrease in 
pressure off axis. The pressure is calculated from the relation 

p - P Q ~ 2 

5 * =(ye , which is obtained by integrating the 0 momentum equation. 

1/2 p V 

* 00 00 


Here P is the pressure on the body surface, P q the stagnation pressure. 
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1/2 p^V^ 2 , the freestream dynamic pressure, and & , the so-called 

pressure relief function which is calculated numerically. 

At the same time Kemp (Ref. 7) presented a similar theory with 
similar assumptions but further restricted his analysis to R em « 0 and 
small interaction parameter. His initial results show a pressure in- 
crease but the analysis was later questioned by Freeman (Ref. 8) , [See 
also Author's Reply (Ref. 9).] The conclusion of this work was that the 
pressure was in fact unaffected for small interaction. 

The first experimental work, due to Zeimer (Ref. 10) investigated the 
alteration to the shock standoff distance with field. This investigation 
was conducted in an electromagnetic shock tube facility using a model 
with a transient magnetic field yielding results over a substantial range 
of interaction. Reasonable agreement for shock standoff distance was 
obtained when compared to Bush's results. 

Wu (Ref. 11) considered constant viscosity, density and conductivity 
model in an analysis similar to Bush. His results show that the effect 
of viscosity is to decrease the magnitude of the alteration to shock 
standoff with field. Lykoudis (Ref. 12) treated a flow model similar to 
that of Bush but assumed the Newtonian- Bus emann pressure relationship in 
his solution. This permits him to obtain an analytical Solution for a 
small range of interaction because the pressure distribution (on the 
body, in particular) is specified. He also assumes a constant radial 
magnetic field which is unrealistic physically, inasmuch as a practical 
magnet contained in the body would have a decreasing radial field. To 
account for this fact the author has suggested a mean value for the field 
in Ref. 13. 
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Pai and Kornowski (Ref. 14) have performed the same analysis as Bush 


except that the magnetic field specified at the shock is assumed to vary 

-3 

as some power of the radius (r for a dipole). The solutions presented 
for different values of the exponent show the same qualitative result with 
regard to the pressure relief as Bush's solution. 

In Ref. 15 Levy has treated a simplified two-dimensional magnetoaero- 
dynamic flow over a cylinder. An attempt is made to account for the Hall 
effect by reducing the scalar component of conductivity in the direction 

a 


of V x B as follows: o 


red 


1 + C 


*2 where is the Hall coefficient. 


H 


This, of course, does not account for the transverse Hall currents. The 
importance of the Hall effect in reducing the interaction is emphasized 
by his simple example. 

Chuskin (Ref. 16) has treated the hypersonic flow around a magnetized 
sphere and cylinder. His model is inviscid but includes compressibility 
and the electrical conductivity is given by a specified temperature de- 
pendent law. The problem is formulated by the method of integral relations 
and then solved by numerical techniques. In formulating the problem in 
terms of integral relations (the method of Belotserkovskii) , Chuskin 
obtains an advantage over the formulation of Bush, in that the analysis 
is not limited to small angle approximations. Variations in physical 
properties behind the shock may be easily introduced. Results obtained 
show that as is the case when viscosity is included, allowing the conduc- 
tivity to vary causes a decrease in effect. His results are shown only 
for relatively small interaction and in the limit of = <», While the 
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standoff Increases with field so does the pressure off axis on the body 
which appears to be contrary to Bush's results. Heat transfer is also 
found to increase. 

Smith, Schwimmer and Wu extended earlier works (Refs. 11 and 18), 
which included viscosity in an analysis similar in method to that of Bush. 
This work was questioned by Porter and Cambel (Ref. 19) on the inclusion 
of viscosity and certain details concerning the evaluation of boundary 
conditions. Correction of these errors resulted in agreement with 
previous works (Ref. 20). 

A group of papers (Refs. 21, 22, 23) deal with the opposite extreme 
from that discussed in previous works. Here the magnetic interaction is 
taken as dominant, the shape and in fact the presence of the body make 
no difference. Levy and Petscheck (Ref. 21) first analyzed the case of 
two-dimensional hypersonic flow over a straight current carrying conductor. 
The magnetic Reynolds number was first assumed small due to small conduc- 
tivity in a preliminary analysis. Subsequently, the high Reynolds 
number — high conductivity case was treated. The density was taken 
equal to pressure as the thermodynamic relationship. The results show 
an increase in shock standoff with increasing wire current. At high 
currents the shock layer was found to lift completely off the wire. 
Experiments (Ref. 22) were performed in a Joule-heated, helium-driven 
gas, diaphram shock tube. The results indicate reasonable agreement with 
the theory for shock position, although the photographs presented do not 
indicate a no-flow region adjacent to the body, a feature of the theoretical 
analysis . 
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A later and more detailed analysis of the corresponding three-dimen- 


sional problem was made by Levy et al. , (Ref. 23) .The existence of a 
magnetically dominated flow for this problem assuming negligibly small 
magnetic Reynolds numbers is shown. The effects of variable properties 
in the shock layer are treated although the Hall effect is not. The 
analysis shows that the shock layer occurs at some position ahead of the 
magnetic source so as to balance the magnetic forces and the momentum 
change of the incoming flow. Behind the shock layer and in the stagna- 
tion region a region of "slow flow" forms which rejoins the shock layer 

outside the stagnation region. These results were obtained in the limit 

1 /^i P 

of e small where e is the density ratio across the shock — . In 

Ps 

considering a problem in which this restriction may not be realistic, it 

is concluded that similar phenomena should result but less distinctly. 

It is also commented that the unexplained differences in the results of 

Refs. 6 and 7 for the pressure terms are resolved by considering the 

location of the source relative to the shock and shock layer center of 

curvature, although these statements appear to be based solely on the 

original papers and not on subsequent work. 

Ericson and Macuilatis (Ref. 3) have applied an analysis similar to 

Bush to the problem of magnetoaerodynamic lift. While primarily concerned 

with the application of this phenomena to flight, several parts of the 

analysis are pertinent to this discussion. In choosing the "volume of 

interaction," the volume of flow which contributes substantially to 

magnetoaerodynamic forces, these authors take 0 = 60° (the angle at 

c 

which the calculations are terminated) to account for the outward move- 
ment of the sonic line when the flow is slowed. Although the Hall effect 
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is not generally treated they utilize a reduced conductivity similar to 
that of Ref. 15 in evaluating the magnetic interaction parameter. In 
later work (Ref. 24) the analysis is reviewed and the magnetically 
dominated flow concept is discussed. 

In an extensive investigation. Porter (Ref. 25) has calculated 
magnetoaerodynamic flow around a hemispherical ly capped cylinder. Using 
a constant property model and like Bush, specifying the magnetic field 
at the shock, he was able to numerically obtain solution for the flow 
field, the shock standoff, the drag coefficient and the pressure distribu- 
tion. The Hall effect, the effect of viscosity and the effect of finite 
magnetic Reynolds numbers are taken into account in separate solutions. 

To obtain the drag coefficients it was necessary to specify a volume of 
interaction outside of which it is assumed that no interaction takes 
place. His results for the shock standoff reduce to those of Bush for 
no Hall effect or viscosity. Likewise, the pressure is found to be de- 
creased for this case, but it is significant to note that for a constant 
value of the interaction parameter the pressure relief is seen to lessen 
and finally become a pressure increase as the Hall effect is increased. 
Although he was not able to handle viscosity and the Hall effect simultan- 
eously, it is suggested that the viscous solution merges with the invscid 
solution for large interaction (S = 50); Thus, an extended composite 
theory is offered to account for these effects together. 

Seeman and Cambel (Ref. 26) have reported an experimental investi- 
gation of total magnetoaerodynamic drag and shock- standoff for a hemis- 
pherically capped cylinder (conducting and nonconducting materials were 
both used) with a self-contained electromagnet. Tests were conducted in 
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a continuous arc-heated facility. Increases in drag and shock- standoff 
distance were measured for interaction parameter (defined at the shock) 
values of the order of S g = 1. These results are presented in terms of 
applicable dimensionless parameters and the standoff measurements are 
compared to the theories of Bush and Lykoudis. Agreement with these 
theories is found by multiplying the magnetic interaction parameter by 
the shock density ratio, e • Drag increases of the order of 10% were 
obtained and a qualitative evaluation of the drag data was made. More 
recently R. Nowak, et al . (Ref. 27) have investigated the same problem 
as Ref. 26, but extending the conditions of the experiment to a wider 
range of pressure, density and interaction. Total drag measurements for 
several flow conditions are brought together with the parameter Se . The 
results show agreement with the theoretical predictions of Ref. 25, thus 
emphasizing the influence of the viscous and Hall effects in reducing 
total drag and shock standoff distances. The results of Refs. 25, 26, 
and 27 have recently been discussed further by Cambel (Ref. 28) along 
with a review of the entire problem of blunt body magnetoaerodynamics. 

1.2 Description of Magnetoaerodynamic Flow 

In reviewing the literature we have seen that the gross qualitative 
aspects of this problem, the change in shock standoff and total drag on 
the body, with increasing magnetic field are reasonably well understood 
in theory and are confirmed by experimental observation. The qualitative 
features of the alteration to the pressure distribution, however, do not 
appear to be in agreement in various analyses and no experimental measure- 
ments of the pressure distribution have been made. 
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We have also seen that most analyses up to the present deal with 
highly idealized models due to mathematical complexities. The limita- 
tions of theoretical treatments may be grouped into three general areas: 
the initial assumptions made; the expected validity of the results; 
and the ability to handle several complicating factors simultaneously. 

A major simplifying assumption commonly made is constant proper- 
ties in the shock layer. While the reasonableness of this assumption 
depends somewhat on the physical example to which it is applied, one 
could immediately question the assumption of constant conductivity as 
Chuskin (Ref. 16) has done. This assumption also requires the concept 
of a "volume of interaction" in order to calculate the force on the 
magnet or the pressure force (Refs. 6 and 25). 

The usefulness of theories which are restricted to small interaction 
is limited to predicting trends and ordering the magnitude of effects. 

Even when the interaction is taken as arbitrary (for example, Ref. 25) 
the validity is often restricted to stagnation point analysis, with the 
accuracy decreasing as the region is extended. Of the combination of 
effects which may be present those of viscosity and Hall effect and 
possibily compressibility have not been treated together, although as we 
have noted earlier, an approximate theory for considering the Hall effect 
and viscosity together has been presented in Ref. 25. Elementary calcula- 
tions have shown that these effects will occur together in the laboratory 
and some flight regimes (Ref. 27). We are therefore limited to the 
knowledge that the inclusion of dissipative mechanisms generally leads 
to a decrease in magnetoaerodynamic effect. 
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To discuss the analytical and physical aspects of the problem we 
require a realistic model, yet one not so complicated as to be unwieldy. 
Bearing in mind that many of the theoretical limitations previously 
discussed are those features which also make the problem manageable, 
we write the following equations governing the problem and make some 
preliminary assumptions: 

Time independent Maxwell equations; 


7 . B * 0 



7 x E ■ 0 


excess charge density, £q = dielectric constant) 

(l.D 


7 x B = ^ J (the permeability of the plasma equals 
that of free space) 


Conservation of Charge: 

V . J = 0 (1.2) 

Ohm t s Law; 

In simple scalar form, J = <j(E + V x B) (1.3) 

Including the Hall effect Ohm's Law may be written as a vector equation 
(Ref. 15), 

J = a(E + V x B) + C (J x B) (1.4) 

H 

where a is the scalar conductivity. An alternate (but not equivalent) 
expression takes g as the conductivity tensor (Ref. 29). Thus, 

J = a • (E + V x B). (1.5) 
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These expressions assume the momentum uncoupling between the ions and 
atoms (ion slip) to be negligible. 

The applicable magne togas dynamic equations may be written as 
follows for steady flow (Ref. 28). 

Continuity: 

V • V - 0 (fdr incompressible flow) (1.6) 

Momentum: 

7 r 2 

p — 2 pV x 7 x V = -Vp + rp v + J X B, (1.7) 

2 

taking V for the divergence of the stress tensor, V • t (Ref. 30). 

It is well known that if the constant property assumption is made 
then it is unnecessary to write an energy equation. For the present we 
shall also make this assumption. However, we recognize that a close 
examination of the flow is likely to necessitate an energy equation and 
while a global energy equation may be written as for continuity and 
momentum, it will prove most practical to use the species energy equation, 
in particular for the electrons, since properties such as conductivity 
and radiation from the flow are governed most strongly by the energy of 
the electrons. The electron energy equation will be considered exten- 
sively in Chapter V. 

We now discuss the scaling of the momentum equation, specifying that 

the magnetic field is that of a dipole 

B 

B *= B cos 0 r' + — ^ s i n q ^ (1.8) 

o z 

and remains undistorted. This simplification is in accordance with 
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Ref. 5, if the nondimens ional magnetic Reynolds number 


Rem «* eu a V <0.1 
“ 00 


(1.9) 


where the product e V represents the velocity in the shock layer. 

00 

To nondimens ional ize Eq. (1.7), the velocity and density are 

2 

referenced to their freestream values p 9 V and the pressure to 1/2 p V . 

r OO 00 r 00 00 

The magnetic field is referenced to B q , the stagnation point value and 
the radius r to the body radius R^. Then Eq. (1.7) becomes (taking 
p*, p*> V* and B* as non-dimensional quantities) 


p* V V* Vp* 

- p*(V* x V x V*) - - -s- + 


p V R - p V 

r oo 00 r 00 00 


J* x B* 


( 1 . 10 ) 


We recognize 

Re (Reynolds Number) = 



S (Interaction Parameter) 


° Vo 2 


( 1 . 11 ) 

( 1 . 12 ) 


which together with the non-dimensional Hall coefficient from Eq. (1.4) 

C H = 6 T ei B o (1.13) 

m 

e 

are the important non-dimensional parameters of the problem. 

The interaction parameter is often interpreted as the ratio of 
2 2 

Lorentz force * a B q LV to dynamic force = pV . The Lorentz force should 
be evaluated in the shock layer however, and the dynamic pressure is 
evaluated in the freestream. Observing that we may divide equation (1.10) 
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by p* 


.CL. 

Poo 


the non-dimensional interaction parameter behind the shock 


Se 


«W C 


V 

00 


since e “ — ■ p* across the shock. We may then interpret as the 
velocity in the shock layer. Thus, the parameter e is also important in 
the problem and Se rather than S is the more natural parameter. 

At this point it is well to review the problem, outline and summarize 
the constraints already made. We have limited the investigation to an 
axisymmetric conducting body containing an electromagnet with both the 
magnet axis and body axis aligned with the flow. A sufficiently low 
magnetic Reynolds number has been assumed so that the magnetic field 
which is taken to be a dipole remains undistorted. This parameter re- 
mains to be evaluated for laboratory conditions. For the present a 
physical model based on constant properties has been adopted since this 
appears to have permitted the most extensive theoretical treatments. 

The validity of this assumption is left open for later review. 
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CHAPTER II 


CONSIDERATIONS IN THE MEASUREMENT OF DRAG AND SHOCK STANDOFF 


The preceding literature survey indicates the lack of experimental 
work pertaining to magnetoaerodynamics in contrast to a number of 
theoretical treatments . The purpose of this chapter is to discuss labora- 
tory investigation of this problem. First, quantities which can be 
measured are formulated from the flow description made in Chapter I. 
Previous techniques are reviewed with emphasis on technique and the 
limitations of these works is used to provide direction for this investiga- 
tion. Finally, preliminary experiments made in this investigation are 
discussed and the shortcomings of these are noted. 

2 . 1 The Components of Magnetoaerodynamic Drag 

The total drag on a blunt body in magnetoaerodynamic flow is composed 
of the pressure force on the frontal surface, the pressure force on the 
base, the viscous force on the body and the reaction force on the magnet. 

In addition to measurements of the total drag, it is desirable to examine 
independently as many of these components as possible, not only for 
individual comparison to theory, but also to understand more completely the 
composition of the total drag. 

Furthermore, the various drags are integrated quantities. If measure- 
ment of the pressure distribution, which is a local quantity, is made then 
additional information concerning magnetoaerodynamic interaction will be 
made available. 

The coefficients of magnetoaerodynamic drag have been formulated by 
Porter (Ref. 25) where the coefficient of a particular component of drag 
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is given as 


Component Force 
D com P ’ 1/2 Pco V co A proj 


( 2 . 1 ) 


where A p ro j 85 » the projected area of the body. We shall use a 

similar development to derive expressions for the component forces. 

Consider the force on the magnet first. This is given as the 
reaction to the component parallel to the freestream of the Lorentz force 
on the fluid and has magnitude 

L f ■ Ivol x | dvo1 <2,2) 


The precise choice of the volume over which this integration takes place 
is a matter of some concern as we have pointed out. To review, theoretical 
treatments, taking a constant property approach and considering the free- 
stream to be nonconducting, take this volume to be a spherical cap on the 
nose of the body, the thickness of the shock standoff. The calculation 

is then terminated at some angle 0 to account for the fact that properties 

c 

such as electrical conductivity are not constant and the influence of 
magnetic interaction is finite in extent. Values of 0 c chosen range from 
0^ = 30° to 0^ = 60°. For the present we will assume 0^ = 45° (following 
Ref. 25). While 0 £ is not a physical quantity, calculation of such quan- 
tities as the force on the magnet are sensitive to the value of 0 c chosen. 
Thus the equation becomes 

R g e c 

L - J J [- (V r B 0 2 - V 0 B r B 0 ) cos 0 + (B r B 0 V r - V^ 2 ) sin 0]r 2 sin 0 drd0 
*** ° (2.3) 
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The viscous drag is proportional to the normal gradient of the flow 
tangential to the body. For the nose region the component of this 
force parallel to the freestream flow is given by 




. 2 
sm 


0 d0 


(2.4) 


where r\^ is the viscosity in the region behind the shock. The shear 
forces on the after body may be similarly formulated. 

Since the pressure distribution is discussed extensively in later 
sections 9 the formulation and theoretical treatment will be examined in 
detail. The 0 component of the momentum Eq. (1.7) is 
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Using suitable functional forms for the velocities and currents the 
authors of Ref. 6 and Ref. 25 find that this equation, which contains the 
pressure, may be directly integrated with respect to 0 under the approxi- 
mation for the stagnation region that 
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COS 2 0 s 1 - 9 2 (2.6) 

sin 2 0 = 0 2 

These particular functional forms are given by (Ref. 25) 

V « 2 e V cos 0 

r oo z 

X 

V fi = -e V Ilisi sin 9 

0 00 X 

V, = e V S i n 0 

ff 00 X 

J, = e V a,B — sin 0 (2.7) 

$ co 1 o X 


where x = ~ , This integration results in 

P - P =| 0 2 @p V 2 (2.8) 

o Z r oo oo 

where & the pressure relief parameter on the body surface (x = 1) is 
given by 

pin 

F'(l) 2 - G(l) + — 

To determine this the pressure distribution on the body, the functions F 1 , 
G, L 1 must be evaluated on the body and suitable boundary conditions 
applied. By application of the governing equations further reduction of 
Eq. (2.9) is possible. We find: 

(inviscid, Rem = 0) ® = e {f' 2 + G 2 - SF* - C^'g} | body 


^ - SL(1) '] 


(2.9) 
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(viscous, Rem =0) ® *= € ^F 1 ^ + F,, '/ Re “ 25F 1 + ^ G u j> |^ od y (2.10) 

Both expressions account for the Hall effect and involve only velocities. 

(P may be determined for the inviscid case (the first of Eq. (2.10)) by 
substituting numerically calculated values of F 1 and G. No such solution 
exists for the second of Eq. (2.10), since this includes both Hall and 
viscous effects. Bush and Porter have determined the pressure relief 
parameter (P for the inviscid case and Porter has extended these results 
to the inviscid case with Hall effect as well as the viscous case without 
Hall effect. The pressure force is then determined by integrating the 
pressure distribution over the nose area of the body. Thus, 

pTT/2 2 

Pf ® 2 J P (0) ttR, sin 0 cos 0 d0 (2.11) 

o 

The magnitude of the base pressure may be calculated using the 
pressure in this region integrated over the area. 

B- = P. A. (2.12) 

f base Dase N 

Since the calculated value of the base pressure is small (Chapter IV) we 
asstime the alteration to it may be neglected and simply correct for the 
no field base pressure in formulating the total drag. The total drag is 
then the sum of these forces: 

Total Drag = V f + P f + B f + L f (2.13) 


2.2 Previous Experimental Work 

Seeman (Ref. 31) has measured the total magnetoaerodynamic drag on a 
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hemisphere capped body in a continuous arc-heated low density facility* 

To accomplish the measurement he used a strain gauge system fixed 
directly to the supporting sting. To provide a magnetic field he used 
a copper wire coil wound on an iron core giving a maximum field of about 
2,000 gauss. This field strength limited the maximum interaction para- 
meter at the shock to about 1 and prevented investigation of the Hall 
effect. Based on theoretical considerations of this problem (Ref. 25) 
we would further conclude that neither the pressure force nor the viscous 
forces were substantially altered when the magnet was turned on so that 
the total drag alteration was due to the force on the magnet. 

Seexnan made preliminary design estimates of this force by assuming 
no Hall effect and using the simple form of Ohm's Law in Eq. (2,3). The 
total force on the body was calculated by a conventional aerodynamic model 
and the heating rate from the high temperature plasma to the body was 
estimated to plan for cooling the model. 

Both Seeman and Zeimer (Ref. 31 and 10) measured the alteration to 
shock standoff distance by direct measurement of photographs. Seeman and 
Cambel have commented on the discrepancies involved in these measurements 
due to the diffuse nature of the shock front (Ref. 26). 

In summary, we expect technical problems similar to those of 
Ref. 31 since our experiment is very similar. For preliminary design 
estimates of drag and heating rates the values given in Ref. 31 will be 
taken. In considering these experiments we would attempt to extend the 
range of interaction and to find a more acceptable method for investigating 
the shock standoff alteration. 
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2.3 Preliminary Experiments 


The arc-heated low density facility constructed for this investiga- 
tion is fully described and analyzed in Chapter III and Appendix A. 

For purposes of the present discussion we note that this is a continuous 
facility as contrasted to typical shock tube facilities which have very 
short duration test times. Thus, the planning of experiments concerning 
magnetoaerodynamic drag need not be limited by the operation of the 
facility but rather by the high heating rates encountered in this type 
of facility as well as limitations on the model itself. In particular 
the magnet test time may be limited by Ohmic heating of the coil. The 
planning of experimental equipment must also take into account the low 
pressure and electrically noisy environment found in these types of 
facilities. 

With the background provided by Ref. 31 a preliminary drag measuring 
system using semiconductor strain gauges was constructed and tested. 

This was found to be unsatisfactory due to hysterisis in the mechanical 
system. Instead a linear differential transformer was found to give 
good results once adequate shielding against the high enthalpy flow had 
been installed. Two complete mechanical systems were designed and con- 
structed, and this final design will be described fully in Chapter IV. 

The only important difference between these two designs is that efforts 
were taken in the second system to keep the sting support small in com- 
parison to the body diameter to minimize interference with the base region. 

In considering Seeman*s experiment it seemed to be most important 
to increase the magnetic field strength in order to increase interaction. 

An obvious advantage is that the magnetic field enters into the interaction 
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parameter as the second power in Eq. (1,12). 

The use of superconducting magnets was considered but these were 
found unsatisfactory after design estimates were made because the thickness of 
the Dewar flask (for cooling) surrounding the magnet lowered the field 
strength to less than that attainable by conventional magnets. Through 
experimentation it was found that the field strengths of conventional 
magnets were increased by the choice of better core material and careful 
positioning of the magnet with respect to the nose of the model. Tests 
also indicated that for the duration of the drag tests (~ 2 sec) the 
wire could withstand a considerable power overload. Final design and 
testing of these magnets is also described in Chapter IV. 

Preliminary drag measurements were conducted using a three inch 
diameter model positioned at 12 1/2" from the nozzel exit. While sub- 
stantial increases in drag could be measured, the flow was unsteady at 
this position giving rise to scatter in the data due to mechanical noise. 

Also it was decided that the model was too large for the flow. Experi- 
ments were conducted to determine whether these problems could be 
corrected by varying the flow conditions either by changing the nozzel 
configuration or varying the test conditions. Although some improvement 
was noted, it was found to be much better to move the model to 4" from 
the nozzel and use a 1 1/2" body diameter. 

Total drag experiments using a 1 1/2 H diameter body at this position 
over the range of flow conditions available in the facility. This initial 
investigation is described in Ref. 27. 

During this work it was found that the flow conditions available in 
the facility were limited either by the facility itself or by the fact 
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that at an available condition the magnetic interaction was too small 
to give substantial results. In summary, this preliminary work served 
to establish the conditions of the flow and facility under which 
experiments could be conducted and indicated directions for the de- 
sign of final equipment. 
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CHAPTER III 


THE EXPERIMENTAL FACILITY AND EVALUATION OF FLOW CONDITIONS 


The purpose of this chapter is to present background material for 
the experiments and discussions in subsequent chapters. The facility has 
already been mentioned briefly and is the subject of the first section of 
this chapter. The facility chosen for this investigation was an arc- 
heated, low-density plasma wind tunnel, similar in nature to the type of 
facilities used by other investigators in experiments requiring high 
enthalpy, hypersonic flows and reasonably long test times (see for 
example, Refs. 32 and 33). 

Analysis of the flow conditions in this facility is required next 
and serves two purposes. Through diagnostic measurement of flow conditions 
the nondimensional parameters appropriate to the experimental problem may 
be evaluated, the most important of these being the interaction parameter 

S, the Hall coefficient C and e, the density ratio across the shock. 

H 

Furthermore, the conditions of the flow where tests are to be conducted 
may be determined. This includes not only measurements of flow uniformity 
but also calculation of parameters such as viscosity and determination of 
the chemical state of the flow. 

The final section of this chapter contains a discussion of the con- 
sistency of the data and calculated quantities and an evaluation of the 
error involved in determining the important nondimensional parameters. 

3.1 The Experimental Facility 

The Gas Dynamics hyperthermal facility consists of a commercial 
direct current arc-heater, which exhausts into a 5 x 11.5 ft. cylindrical 
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vacuum tank. The expansion of the gas is through a conical nozzle 
1 1/2" diameter exit with a 9/32" diameter throat and a half-angle of 
15°. 

The vacuum tank is continuously maintained at a low pressure by a 
5,000 cfm impeller blower which is backed by a 460 cfm rotary vacuum pump. 
The pressure is read by an "Alphatron" Gauge. Inputs for water and 
electricity are through the sides of the tank. The entire facility is 
controlled from an operating console and the flow is visible from this 
station through a viewing port. The facility is shown pictorially in 
Figure 2 and photographically in Figure 3. A more complete description of 
the facility is contained in Appendix A, along with specifications for 
the various equipment mentioned here. 

Auxiliary equipment for use in this facility includes an overhead 
table which, controlled from the console, can be moved in three dimensions. 
This is used to support both the diagnostic devices and the models under 
test. Position of the table is read by cam operated microswitches and 
indicated at the console. Also available for the facility is a fast-acting 
heat shield which when in position deflects the flow to the side and allows 
a model or probe to cool. This is also controlled externally. Power for 
the arc is taken from the common laboratory power supply consisting of four 
rectifiers, commercial welders. The output may be filtered if required. 

The external variables were mass flow and power which is measured by 
reading the voltage across torch and current through it. Heating losses 
were determined by measuring the flow rate of the cooling water and the 
temperature rise by thermocouples. The pressure in the stagnation chamber 
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and arc chambers were variously read with mercury manometers and commer- 
cial pressure gauges* 

3.2 Choice of Flow Condition 

For experimentation we require a uniform continuum flow in which it 
is possible to produce high interaction. Data on the operation of the 
facility is presented in Appendix A. The operation of this facility is 
restricted to what may be characterized as a viscous to merged layer 
regime due to the available pumping capacity. Furthermore, it is diffi- 
cult to change any of the external variables of the facility without 
substantially affecting all the flow quantities. (Cf. Figure A2 which 
shows some operating parameters of the facility plotted against mass flow 
for a constant arc current. For example it is seen that the mass flow 
cannot be changed without simultaneously changing the ambient pressure 
and the power input.) Over the limited range of variation where the 
flow is continuum, viscosity is always important. 

We now consider the flow uniformity and interaction. The maximum size 
of the model is dictated by the maximum flow size. The model may be made 
smaller but this decreases the interaction parameter not only since R^ 
decreases but also because the maximum stagnation field strength de- 

2 

creases. This is more important since the interaction depends on B q . 

The model size selected on these considerations was R^ = 3/4". Then the 
conditions which seemed to give the most uniform, high conductivity flow 
were selected as the single operating conditions. (Some samples of radial 
pressure profiles at other conditions are shown in Figure B-2. For example, 
it is seen that moving the axial position of the pressure probe while 
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holding the flow rate constant results in a more nonuniform flow. 

Changing the flow rate at the running position also increases the non- 
uniformity,) In addition, only argon gas was used, being attractive due 
to the low ionization potential, chemical inertness and large amount of 
data available for this gas. 

It is realized that chosing a single flow condition is a severe 
limitation on the experiment. However, we emphasize that if different 
running conditions and model sizes were chosen new problems in relating 
results would be introduced. As was mentioned above, since no substantial 
change in such important parameters as Reynolds number can be made, no 
advantage is gained. 

The running condition chosen for this investigation is most conven- 
iently specified by stating the external quantities. These are: 

^Torch = 960 am P s 



Under these conditions the supersonic stream is self-luminous and Mach dia- 
monds may clearly be seen. The flow is over expanded at the test station. 


28 



3.3 Diagnostic Measurements of Flow Conditions 


The research program of which this investigation is part includes a 
series of diagnostic measurements of flow conditions. While a number of 
techniques are available as a result of current research, we have concen- 
trated on two probes of the type used in gas dynamics flow research, 
and spectroscopy. Other techniques were found to be either inconvenient 
because of the facility design or could not be used because of the flow 
conditions. Spectroscopic measurements of electron temperature and 
number density will be fully discussed in Chapter V. For the present we 
will be concerned with the probe measurements of pressure and mass flux. 
Additional information and specifications may be found in Appendix B. 

For the purposes of pressure measurements, a differential pressure 
transducer was mounted in a water-cooled chamber suspended from the over- 
head table. From this chamber a water-cooled sting extends down into the 
flow and different shaped heads may be soldered on the sting. The 
assembly is shown in Figure 4. This apparatus was used for a variety of 
impact, wedge, and static pressure measurements. 

For impact pressure measurement, blunted cones of 15° half-angle and 
a hemispherical shape were used with 1/16" and l/8 n taps. Static pressure 
measurements were made with a 1/2" diameter probe with taps normal to the 
surface located 10 diameters from the tip. Both 1/8*' and 1/4 M taps were 
used. 

A mass flux probe was designed and constructed. The principle of this 
probe is that instead of forming a bow shock ahead of the probe the shock 
is swallowed and stands inside the probe. The fluid from the freestream 
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which would normally flow around the probe is trapped in the probe. Thus, 
the mass flux into the probe represents the freestream mass flux. In 
order to accomplish this the probe was constructed in the following 
manner. A sharp-edged opening (both 1/8" and 1/4” diameter were used) is 
located at the tip of a 15° cone. The diameter of the inside channel then 
enlarges and then exhausts into a large reservoir vacuum tank. This 
reservoir is maintained at a pressure lower than the tank pressure before 
the test. The flow into the reservoir is initiated by opening a control 
valve. After a fixed period of time the pressure and temperature are 
measured at the reservoir, thus the mass flux can be deduced for choked 
flow. 

The results of probe diagnostic measurements are as follows: 


^impact 
P „ 

static 


6.4 mm 

20 ^ 


P v 

'00 00 


= 8.7 x 10" 3 gm/cm 2 


sec 


(3.2) 


Of these, the impact pressure appears to be most reliable. The static 
pressure is definitely lower than the ambient pressure of 400 but the 
low Reynolds number makes the measurements difficult to interpret. The 
measurements are also in question due to the difficulty of insuring 
that the shock is truly swallowed (choked flow) and the question of the 
influence of viscosity on the measurements. A detailed discussion of 
error is deferred until later, but we note here that the devices such as 
the pressure transducer typically have an error of the order of one per- 
cent. It would be optimistic then to take this as the error of the 
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measurement since scatter in the readings and system error is larger. 

From the results of Chapter V we have the electron density in the 

l^h O 

freestream as 1,7 x 10 /cc and the electron temperature as 5,000 K. 

Sample radial profiles of n^, and impact pressure are shown in Figure 5, 
At this point we may discuss the choice of model test station. For 
the running condition chosen the model was positioned approximately in 
the center of the first shock diamond (4 n from nozzle exit) which was 
found to be the most uniform area both axially and radially. 

3 .4 Calculation of Flow Quantities 

Using the measurements just discussed and the spectroscopic results 
reported in Chapter V we proceed to calculate the various flow quantities 
which will enter into discussions of the results of this investigation. 

This analysis will be guided by the results of recent investigations by 
other workers using low density arc-heated facilities (see, for example, 

Ref. 32). A well known characteristic of this type of facility is that 
the freestream is out of equilibrium in the sense that the ions and atoms 
have a lower kinetic temperature than the electrons. This is due to the 
fact that the electrons and heavy particles are not well coupled by energy 
transfer between species. The high electron temperature tends to be 
maintained during expansion while the heavy particle temperature drops in 
the supersonic freestream as in normal supersonic flow. 

It is reasonable to assume however, that equilibrium is approached in 
the arc and stagnation chamber due to the relatively high pressure (0.5 atm) 
and slow flow velocities. From the net input power of 7.7 KW to the gas 
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and the flow rate of 1.23 g/sec we may calculate the stagnation enthalpy 
to be h Q « 2,700 BTU/lbm. From the equilibrium calculations of Arave 
and Huseley (Ref. 34), and knowing the stagnation pressure we find for 
equilibrium argon in the stagnation chamber 


T - 9, 700°K 
o 

a - .025 


(3.4) 


where Ot is the degree of ionization. 

With knowledge of the electron temperature in the freestream we 
assume that the flow is chemically frozen (i.e., that recombination in 
the freestream is a minor effect and since the temperatures are low 
subsequent ionization is not important) . This means that the ratio of 
specific heats, y, is simply that of a monatomic gas, 1.67. On this 
basis using the ratio of impact pressure to static pressure (= 33.5) the 
tables of Ref. 35 show for isentropic flow 


for which 


hence 


n. - 4 - 75 

1/e - ~ - 3.5 

Pm 


T = 1100°K 
00 


(3.5) 

(3.6) 

(3.7) 

(3.8) 


where T is the heavy particle freestream pressure. 

OO 

The number density of the freestream is more difficult owing to 
interpretation of the n V measurement. Instead of using this measurement 

“oo 00 
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directly, we shall estimate a set of flow-averaged quantities in the 
following manner, for a freestream heavy particle temperature of 1100°K 
the equilibrium data (Ref. 34) for the speed of sound gives 


a ** 2100 ft/sec = 6.4 x 10 cm/sec 
00 


from this we conclude 


v = M a = 3 x 10 J cm/sec 
00 00 00 


(3.9) 


(3.10) 


Taking an appropriate jet radius from probe and spectroscopic data (see 
Figure 5) we may estimate the area of the jet to be 


A. . = 37 cm 
jet 


(3.11) 


Therefore o V = m/A 

“no m 


jet 


-2 2 

3.35 x 10 gm/sec cm . Hence from our previous 


estimate of v 


-7 


p = 1.1 x 10 g/cc 

r OO 


(3.12) 


From this an average value of the electron number density in the free- 
stream may be calculated 

n = a — = 4 x 10 13 /cc (3.13) 

eo° m 

h 


where m is the mass of an argon atom. In order to make a better evalua- 
h 

tion of o V it is now possible to present another method of obtaining 
r 00 00 

this quantity. From the results of Chapter IV we may introduce two more 
independent observations of p^V^. Values of the total aerodynamic drag 
and pressure drag coefficients with no magnetic field are given by Ref. 25 
as 1.5 and 1.0 respectively. Using the measured values of these drag 
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(3.14) 


forces and the definitions of the drag coefficients 


C drag 


drag force 
1 2 2 
2 Poo V co *b 


the value of p V may readily be calculated. Then using the value of 
r 00 00 

V = 3 x 10 5 cm/sec calculated previously p V is obtained. The values 

OO '00 oo 

calculated in this manner are 


-2 2 
p V = 2.5 x 10 g/sec-cm 

“oo oo 


-2 2 
= 2.7 x 10 g/sec-cm 


(from total drag = 65.6 g) 

(3.15) 

(from pressure drag « 48 g) 


We note that these calculations are also averaged values (although 

the resolution is somewhat improved) since both the drag forces and V 

00 

are averaged quantities. Also obtaining radial variations is not 
possible by this method. The values obtained by this method lie between 
the probe measurements and the flow averaged mass flux. We would expect 
them to lie closer to the latter quantity. 

We now proceed to calculate additional quantities which will be 
needed in later sections. In the following estimates we take 
N = 5 x 10^/cc; N = 10^/cc, T « 10^°K T = 5 x 10^ °K. First 

es > eoo 9 es eoo 

we find the conductivity and viscosity from the equilibrium calculation 
of Ref. 31. 


a , n = 2.5 x 10 mho/m 
shock 

-4 

u - 5.4 x 10 g/cm-sec 

oo 

n shock = 8 X 10 ’ 4 s/ cm - sec 


(3.15) 
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II II III 


For future reference and comparison we now calculate some microscopic 
quantities: The root mean square velocity 


3kT 1 ■ 5 x 10 7 cm/sec (freestream) 

(rzr^) ' 


m 


■ 6 x 10 cm/sec (shock) 


(3.16) 


3kT, 

"h 


7 x 10 cm/sec (freestream) 
2.5 x 10^ cm/sec (shock) 


where the subscript h denotes the heavy particles. The mean free path for 
the heavy particles (Ref. 36, 37) 


Neutral - Neutral 


X 


aa 


16 r m A g/2 t 
5 ' > 2rrkT h '' P 


« 0,13 cm (freestream) 
^ 0.03 cm (shock) 


(3.17) 


Neutral - Ion 


X 


ia 


0.32 


+ 142 

T + 11.5 T h 


aa 


= 0.035 cm 


(3.18) 


The electron-atom mean free path is taken from Ref. 36) to be 


X « 0.2 cm 
ea 


(3.19) 


The Debye length is calculated (Ref. 37) 

x = (-^) 1/2 

*D V 4TT n ' 
e 

J33 10 cm 


(3.20) 
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for the range of n and T 
e» ( 

by Ref. 37 


0) 


c 



estimated. The cyclotron frequency is given 
■ 1.75 x 10 11 B/sec (3.21) 


when B is in Webers/m . Two collision times will be needed: the 

collision time between electrons and ions (Ref. 37) 

T ei » 1.45 x 10“ 10 sec S (3.22) 

and the collision frequency between ions and neutrals 

V. r 

V. ® t — = 2 x 10 /sec (3.23) 

la *ia 


The ambipolar diffusion coefficient is given by (Ref. 38) 

k(T i + T e>-. 4 2, 

D = = 7.5 x lCr cm /sec 

ao ° m i v ia 


(3.24) 


3.5 Calculation of Dimensionless Parameters 

On the basis of the preceding calculations the appropriate dimension- 
less parameters may be evaluated. 

1. The free stream Knudsen number based on the radius of the body 


“b 


0.75 in. 


Kn 


00 



(3.25) 


2. The Reynolds number at the body 


Re 


VH> 

^8 


*5 0 ( 100 ) 


(3.26) 
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3. 


The Mach number and density ratio as previously calculated 


M =4.75 
00 

e E 0.29 


4. The freestream ambipolar Schmidt number 

^00 


Sc 


0 D 

r <x> aco 


0.6 


(3.27) 


(3.28) 


5. The magnetic Reynolds number 


Rem — cV^R^ ^ 0(10 ) 


(3.29) 


The interaction parameter using p calculated from the drag 
coefficients 

2 


2 

a B R, 

S = „ - = 2.0 B 

p V o 

~ m oo 


(3.30) 


(B q the magnetic field at the stagnation point of the body 
expressed in kilogauss squared) 


The Hall Coefficient 
eB 


C„ *= T . = 0.8 B (B in Kilogauss) (3.31) 

II m ex o o 

e 


3.6 Consistency of Data and Error 

The underlying assumption in the analysis of this type of facility is 
that chemical reactions are not of major importance so that the concen- 
trations of the species are relatively constant. This will not be the 
case over large axial distances since the various species will tend to 
equibrate both energetically and with respect to concentrations. However, 
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in the regions of interest, the shock layer and the frees tr earn flow from 

the nozzle, the assumption of frozen flow may be justified since the scale 

lengths involved are short and the flow velocities are high. We return 

to a discussion of the shock layer in Chapter V. The freestream is 

considered now since the assumption of no chemical reactions here lead to 

the value of y *= 1.67, the specific heat ratio for a monatomic gas. 

We expect three-body recombination to be the dominant chemical 

reaction in the freestream. Picking a value of a the rate coefficient 

-10 3 

from Ref. 39 of 1.1 x 10 cm /sec, we have the rate equation 


dn 

e 

dt 



which may be rewritten as 


(3.32) 
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taking Az as 1 cm we have 
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e = io " 2 


(3.34) 


as an estimate of the recombination effects. Since this is the principle 
chemical reaction we conclude that taking y = 1.67 is reasonable. 

The usefulness of y was in the evaluation of the Mach number and 
e, which was also questioned because of the difficulty encountered in 
measuring the static pressure accurately. An independent and direct 
observation of e is available from the results of Chapter V where the 
shock density profile is measured. The calculated value of e = .29 is in 
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good agreement with these results. 

Careful consideration must be given to the value of n V which enters 

•00 CO 

into the interaction parameter calculation. Four independent measurements 
were given, three showing consistency and the probe data being less 
reliable. To summarize the results were 

-2 2 
( ca ^ cu ^- atec ^) = 3.35 x 10 g/sec cm 

o 

Poo V w (pressure coefficient) = 2.7 x 10~ (3.35) 

Poo Vo ( drag coefficient) = 2.5 x 10“^ 
p V (probe) = 8.7 x 10*^ 

'CO 00 

As a final consistency check we observe that the centerline number 

density measured by spectroscopy (1.8 x 10^/cc) may be compared to that 

13 

value obtained by flow averaged calculations (4 x 10 /cc). This agree- 

ment is reasonable allowing for gradients in the flow and the relatively 
low accuracy of the spectroscopic method in obtaining absolute number 
densities . 

In evaluating the nondimens ional parameters in the last section, 
values for the various flow quantities which seemed most accurate and 
consistent were selected. Obviously, other evaluations are possible. In 
particular, it was found most consistent to use the value of n V calculated 

v oo oo 

from the drag and pressure coefficients in the magnetic interaction para- 
meter S, due to problems in interpreting the mass flux probe measurements. 
While an absolute error estimate for S (and the other parameters) cannot 
be made, it is possible to examine the error in S calculated in this 
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manner by estimating the error in the values of used. We assume 

that the error in reading the magnetic field is negligible in comparison 
to these measurements. We take the difference between o V calculated 

r OO CO 

from the pressure and drag coefficients as an estimate of the error in 

P V * 

r 00 00 

We will be primarily interested in the interaction parameter and the 
Hall coefficient in later chapters. Thus, 



where t and n are calculated from the centerline variations of n and 
e e 

T fi . The error in the measurement of magnetic field will be neglected. 
The commonly used difference method for determining the error in a 
calculated quantity gives 


fs 



15% 
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(3.37) 


40 



CHAPTER IV 


EXPERIMENTAL OBSERVATIONS OF MAGNETOAERODYNAMIC DRAG 


It was established- in the first two chapters that the total magneto- 
aerodynamic drag was the sum of several components, each altered by the 
interaction of the flow with the magnetic field. 

Preliminary measurements of total drag alteration such as those 
reported in Chapter II indicated that measurement of the individual 
components of drag would provide better understanding of the measure- 
ments of total drag. The following approach was taken. The total drag 
with and without magnetic field and the force on the magnet were 
measured. Then a detailed study of the pressure distribution on the body 
with and without magnetic field was undertaken, leading to calculation of 
the pressure force. The viscous force, being the smallest component and 
most difficult to measure could then be calculated from these measure- 
ments assuming base drag was unaltered. Thus the total drag and individual 
components were available experimentally for discussion and comparison to 
theory . 

These experimental measurements are presented in the following 
sections along with a description of the apparatus in each. Typical data 
are presented as functions of magnetic field squared, representing a dim- 
ensional form of the magnetic interaction parameter. This is done for 

convenience, since an evaluation in terms of an alternatively defined 

2 

magnetic interaction parameter may be made by simply multiplying B by a 
constant factor. The data will be fitted to empirical relations and 
replotted in terms of the interaction parameter in Chapter VI for purposes 
of discussion and comparison to theory. 
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4.1 Design, Construction and Calibration of Experimental Apparatus 

The bodies and magnets used in these experiments were made as uni- 
form as possible to eliminate systematic error. The basic body design 
was chosen to be 1 1/2" in diameter with the models containing magnets 
having a 1 1/2" afterbody and the pressure models without magnets, 1/2" 
after body. The magnets are fitted close to the front of the body to 
maximize the field as shown in Figure 6. Figure 7 is a photograph of 
a magnet used and a total drag model along with pressure, Lorentz drag 
and flat-end bodies for corapairson. 

The final magnets were based on experimentation with a number of 
jig wound magnets. The rough dimensions were as follows: 1” long by 

1 1/4" diameter with a 3/8 M core. The nose was rounded to insure a 
close fit at the nose, thus maximizing the field. The wire chosen for 
winding was #19 Anaconda HML-coated wire, desirable for its high tempera- 
ture coating. This gauge was heavy enough to take high currents and 
provide good heat transfer when cooling was used yet light enough to wind 
easily. The core material was vanadium permandur, an alloy commonly used 
for this purpose, heat treated according to the manufacturer’s directions. 
The magnets once wound were fitted with cores mounted directly to the 
bodies at the base (see Figure 6). Typical magnet resistance was 0.5 Q. 
Both commercial storage batteries (32) and a commercial power supply were 
used to power the magnets. The batteries provided practically any combina- 
tion of voltage and current desired and produced a very constant power 
for short runs. The commercial supply was used for longer runs, testing 
and charging the batteries. 
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The measurement of field strengths was performed with a Bell 
Gaussmeter, Model 120. The probe area of this instrument is approxi- 
mately 1/8" x 1/4” and of negligible thickness. This provides reason- 
able spatial resolution although not adequate for detailed measurements 
as a function of angle in the stagnation region. 

Two tests were made. The magnetic field strengths at the stagna- 
tion point for various coil currents and the field drop-off with axial 
positions for a constant current were measured. Samples of field 
strength tests are shown in Figure 10 for several models. The normal- 
ized drop-off results are compared in Figure 11. These are compared to 

R -3 

the theoretical dipole given by B = B (— ) and also for an exponent 

° \ 

of 3.6. The dipole field has generally been taken in theoretical treat- 
ments. It is seen that in the shock region the agreement is reasonable. 

To insure that positioning was not affecting the distribution severely, 
one magnet was moved back 0.1” with respect to the stagnation point and 
these results are seen to agree with the previous results. 

To measure the force on the body the following system was constructed 
(shown in Figure 8). The model was mounted on a 1/2” sting carrying 
power for the magnet and cooling water. This sting was suspended by two 
vertical leaf springs so that the force on the model is seen as a small 
displacement of the sting. The sensitivity can be increased or de- 
creased by changing the thickness of the leaf spring. 

This displacement was read by a Sanborn Linear Differential force 
transducer which also contributes to the restoring force on the sting, 
and provides a millivolt D.C. output linear with force. The electrical 
and water input to the sting were arranged so as not to interfere with 
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free movement of the model. Surrounding the whole mechanical apparatus 
is a water cooled shield for protection of the assembly and transducer. 
Further shielding eliminates the drag of the sting and suspension system. 

The system as a whole mounted to the overhead table assembly and the 
model was then positioned by a mechanical jig and a level. Once installa- 
tion was completed, calibration could easily be performed during any 
series of runs. This was done by positioning a small pulley behind the 
model and attaching a thread from the model to the pulley along the 
direction of displacement. The thread then runs over the pulley and 
desired calibration weights attached. The transducer proved to give 
extremely linear response. A sample calibration curve is shown in 
Figure 9. Typical sensitivity was 1.6 MV/g, considered to be adequate 
for these measurements. 

Another important characteristic is the time response of the 
mechanical system due to the short lengths of time during which the 
magnet could be used. Coupled with the time response of the mechanical 
system is the response of the electrical system. The total effect was 
determined by placing a weight about equal to the aerodynamic drag on the 
calibration system. A weight approximately the same as the magnetic 
field induced alteration was then suddenly added and the response measured, 
thus simulating the condition of an actual test. From the time the 
magnetic field was turned on a time of about 0,5 sec elapsed before a 
reliable measurement could be made. This is due partially to the 
mechanical response of the system to a step input and partially 
to the response of the recorder (.5 sec, full scale). With an electrical 
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filter added to reduce the noise the period elapsed before a reliable 
measurement could be made was increased to approximately 2 seconds. 

For the Lorentz forces a similar apparatus was constructed except 
that the body was attached rigidly to the outer heat shield. Thus the 
magnet alone was suspended from the sting and the forces on the magnet 
were measured directly. This apparatus was not water cooled and the 
movable heat shield was employed in front of the body. The calibration 
procedure was the same as that used for the total drag model. 

In order to measure pressure distributions without magnetic field 
solid copper bodies were used, 1 1/2” in diameter with 1/2" afterbodies. 
These were soldered to the impact pressure device described in Chapter 
II. One body was made for each angle measured and 4 taps spaced symmet- 
rically were cut to improve response time. The same pressure transducer 
as in the impact pressure was used. Angles measured in this manner were 
16°, 26°, 32°, 48°, 56°, 90° and the base pressure was measured at 
r = 3/8”. In separate bodies one tap was drilled at 0° and two at 8°. 
Both 1/16" and 1/8" taps were used to determine transpiration effects. 
These models were then suspended from the overhead table and the heat 
shield was used to prevent excessive heating. 

The measurement of pressure distributions with field was more 
difficult due to the lack of space in the body. To overcome this problem 
in the stagnation region, the magnet core was extended through the body 
shell and machined to form the nose. Two 1/16" taps were then drilled 
through the core, one at 0° and another at 8°. 

For off-axis measurements a small copper tube was imbedded in the 
shell and machined flush with the body surface. A tap was then drilled 
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into the tube at the desired angle. Angles chosen were 20°, 32° and 
44°. Only one tap was provided at each angle and several bodies were 
constructed in this manner. All the bodies for the magnetic field 
pressure experiments were attached to a 1/2" sting supported from a 
simple transducer housing. This was then supported from the overhead 
table. 


4.2 Experimental Measurement of Drag and Drag Components 

A typical trace of a total drag test (unfiltered output) is shown 
in Figure 12. The apparent difference in the two traces is due to the 
mechanical crossover in the two channel recorder. The rise of the two 
traces is actually simultaneous. The vibration seen is attributed to 
mechanical noise from the roughing pump. The aerodynamic drag was 
measured at the beginning of a series of tests to minimize errors due to 
heating and these results were averaged, giving 66 g. The total drag 
with field was measured in many different test series over several months 
using different equipment. The results of these tests and the Lorentz 
drag experiments are shown as a plot of percent alteration versus magnetic 
field squared in Figure 13. 

To determine the effects of body shape on the alteration to total 
drag, a similar experiment was conducted with a flat ended body and 
magnet. These results are shown in Figure 14. The results for the 
hemispherically shaped body are also shown in the same figure for com- 
parison. It is seen that the flat ended body shows a much smaller effect. 
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This is not surprising since the aerodynamic drag is higher for the 
flat ended body. The component of drag with field is also smaller which 
is attributed to the fact that the flow is both slower and less perpen- 
dicular to the magnetic field lines around the stagnation region thus 
reducing the Lorentz force. 

The volume of the shock layer however, is greater for the flat end 

body. The influence of this increase as well as the influence of the 

conductivity of the freestream on all drag measurements is resolved in 

the following manner. It is shown in Chapter V that the shock standoff 

for the hemisphere body is about .35”. One may readily estimate that at 
2 

this position B has dropped about an order of magnitude. Although the 
conductivity does not drop much across the shock since the electron tem- 
perature is nearly constant, the velocity vector is axial in the free- 
stream and not perpendicular with the field lines. On this basis it is 
reasonable to assume that neither the freestream nor a slight increase in 
shock standoff contribute significantly to the observed increase in drag. 

The distribution of pressure on a spherical surface as given by 
modified Newtonian theory is 


where 


C = C cos 0 
o 


C P = 


p - p 

*o r 
1/2 p V 

“oo nr 


(4.1) 

(4.2) 


where 0 is the angle measured from the centerline, and the subscript o 
refers to the stagnation point. Accordingly, the first step in studying 
pressure forces is to measure the distribution of pressure without a 
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magnetic field. This will give finally, not only the pressure force 
on the body, but also indirect indications of the effects of flow non- 
uniformity. 

Once the pressure has been measured as a function of angle, the 
total pressure on the frontal surface may be determined by integrating the 
pressure times the projected area upon which it acts. 

The projected area of the frontal surface is given by 

2 

dAproj * cos 0 dA = 2 cos 0 sin 0d0 (4.3) 

so that the elemental force on the frontal surface is given by 

dF = 2P(0) I^ 2 sin 0 cos 0d0 (4.4) 

From this we must subtract the force due to the base pressure to get the 
net force due to fluid pressure on the body. If the base pressure is 
taken as averaged over the area not covered by the sting, this is given by 
the product of the pressure and the area. 

With no magnetic field the pressure distribution was measured with 

both 1/16" and 1/8” taps. The results are shown in Figure 15 normalized 

2 

to the stagnation point value and compared to cos Q . The pressure at 

the shoulder and the base pressure were slightly negative which tends to 

confirm the static pressure value slightly below the tank pressure. The 

1/8” tap gave consistently higher values than the 1/16" tap and also 

2 

followed closely the cos 0 curve. 

The measurement of pressure distributions with field was conducted 
as follows. The pressure alteration was measured as a function of field 
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and angle as shown in Figure 16 for the 1/16" tap size. This data was 
then averaged and added to the no-field distribution for 1/16" tap. It 
was then assumed that the discrepancy between the 1/16" and 1/8" holes 
was due to thermal transpiration and that to correct for this the 
following procedure was followed: we first assume that the ratio of 

pressure with and without magnetic field for the same tap size is 
independent of tap size. Hence, 


P(0 


9 I /O 11 

p ( e »b) 1/16 .. p( 9 ,o) 1/16 „ 


(4.5) 


We now correct for the thermal transpiration effects on a 1/8" tap using 

the results of Ref. 40. Since the correction is small there is negligible 

error in simply assuming a constant value of T^/T^, the ratio of the gas 

temperature at the tap to the gas temperature at the transducer. For 

T, /T =4 this amounts to a 5 % correction, 
h c 

These corrections were applied together to the 1/16" distributions 
and the results integrated giving 48 g the corrected distribution. 

Several of these distributions are shown in Figure 17 (connected to the 
1/8" tap size as shown above) with the most surprising result being the 
pressure rise in the stagnation region. The integrated effect however, 
is still a decrease in pressure force with field as shown in Figure 13. 

The total drag value of 66 g and the pressure drag of 48 g with an 
estimated base pressure of 3 g permits calculation of the viscous force, 
found to be 21 g for the non-magnetic case. Assuming the base pressure 
does not change with magnetic field, the above procedure would allow us 
to calculate the change in viscous drag with magnetic field. 
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CHAPTER V 


INVESTIGATION OF THE SHOCK LAYER 

The shock detachment distance for a blunt body in supersonic flow is 
a quantity frequently measured for comparison to theory . The commonly 
given reason is that this quantity is conveniently measured and a 
theoretical treatment of the blunt body problem must predict this quantity 
If the entire analysis is correct. While few of the theoretical treat- 
ments reviewed in Chapter I consider drag, all determine the shock 
detachment distance in magnetoaerodynamic flows. In theories similar to 
that of Bush (Ref. 6) this quantity is predicted automatically when the 
body has been located during the numerical integrations. 

The visually observed shock layer with and without field may be seen 
in the photograph of Figure 18. Although shock appears diffuse and thick, 
a change in shock standoff is apparent. In this investigation shock 
detachment distance was measured by photography and then a detailed 
examination of the shock layer was made by spectroscopy. The purpose of 
the spectroscopy was twofold: to aid in the interpretation of photographic 

measurements and to provide diagnostic measurements of electron density 
and temperature. Both of these techniques are optical, and so have the 
advantage of not interfering with the field while capitalizing on the 
self-luminous property of the flow. 

5 . 1 Shocks in a Partially Ionized Gas 

The formation of a shock in a partially ionized gas is complicated 
by the presence of three species; ions, electrons, and neutrals, which 
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may be reacting. In addition, nonequilibrium with respect to concentra- 
tion or species temperature may exist and In principle, each species may 
pass through an individual density rise. Several works have attempted to 
treat these problems separately. To provide background for this chapter 
we will review these here. We may then determine the state of the shock 
encountered in this experiment by using the results of Chapter III. 

Talbot and Grewal (Ref. 33) have considered extensively the formation 
of a shock wave in a partially ionized gas where the electrons are at a 
different temperature than the heavy particles. Their approach is through 
the electron energy equation (to be discussed) and assuming no 

energy transfer between the electrons and heavy particles and no reactions 
(frozen flow) . The method of solution is to assume a heavy particle 
shock given by a step function and numerically integrate from the free- 
stream and the region behind the shock and match conditions at the step. 
The fact that the electrons are energetically uncoupled from the heavy 
particles results in a broad zone of elevated electron temperature up- 
stream of the shock forced by the electron compression at the shock and 
transmitted upstream in the electron gas due to the high thermal conduc- 
tivity of the electrons. Electron compression at the shock is found due 
to charge neutrality restrictions which prohibit a change in ion density 
(in this case heavy particle density) without an equal change in electron 
density. Because of the high electron temperature ahead of the shock the 
electron compression is essentially isothermal. These results are offered 
to explain the characteristic "dark space 11 ahead of blunt body shocks 
in low density arc-heated wind tunnels (see Figure 18). The reason 
given for the occurrence of the dark space is that the three— body 
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recombination which determines the radiation from the plasma is decreased 
in this region due to the elevated electron temperature ahead of the 
shock. 


These results were later confirmed in a much more general treatment 
of the problem by Jaffrin (Ref. 41) . He treats the shock structure for 
varying degrees of ionization, from a very weakly ionized gas to a quasi- 
fully ionized gas, by considering the species momentum and energy equa- 
tions as well as the global conservation equations. 

Chen (Ref. 42) has investigated these results in a facility similar 
to the one used in this work, obtaining data similar to the predicted effect, 
i.e., the electron temperature rises before the shock and remains relatively 
constant through the zone of electron compression. A constant electron 
temperature across the shock was also found by Sonnin (Ref. 43) conducting 


an experiment at low gas temperature and ionization (T 


10 


11 


gas 


= 200°K, 


N =— — ). Christiansen (Ref. 44) working with a cesium-seeded argon gas 
e cc 12 

10 x 

heater facility (N fi - 0 ( ~ — ) ) has made a much more detailed study 
involving an investigation of thermal conductivity. Both spectroscopic 
studies and Langmuir probes were used. His results show a small electron 
temperature rise across the shock, so that the temperature ratio across 
the shock cannot be readily predicted. In addition to the analysis of 
charged particle separation occurring at the shock made by Talbot and 
Grewal, Pipkin (Ref. 45) has considered this effect by treating a highly 
idealized one dimensional flow problem. His conclusion substantiates 
that of Ref. 33, as he finds that the electrons tend to diffuse upstream 
so that a charge layer is formed at the shock. The characteristic decay 
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length of this layer is found to be the Debye length. 

Sonnin (Ref. 43) has investigated the question of possible separation 
between the charged and neutral particles by ambipolar diffusion. In 
order to determine the magnitude of this effect it is necessary to 
examine the conservation equation for both electrons and ions so as to 
include the effect of momentum transfer between the heavy particles. It 
is found that when the ambipolar Scmidt number defined in Chapter III 
becomes small the electrons and ions pass through an independent density 
rise upstream of the neutral shock. This result was then confirmed by 
experiment. 

To summarize, the major effects which may be present at a shock in 
a partially ionized gas are: an elevated electron temperature ahead 

of the heavy particle shock due to high electron thermal conductivity, 
separation of the electron density rise from the heavy particle shock governed 
by the Debye length and charged particle neutral particle shock separa- 
tion due to ambipolar diffusion. 

5.2 The Chara.cter of the Shock Layer 

Although phenomena discussed in the preceding section may considerably 
complicate the phenomena in the shock region, a hydrodynamic shock in the 
neutral species is present and we may predict certain features before 
proceeding to analyze the phenomena discussed above. 

Probstein and Kemp (Ref. 46) have categorized the flow through the 
shock by comparing the ratio of mean free path in the shock to body radius 
(shock Knudsen number) . Thus they find that when the Reynolds number is 
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low, 

v*b « * in 
>-j\ « 1 

>-„ /R b < 1 

Using the results of Chapter III, we have 


Viscous Layer Regime 

Incipient Merged Layer 
Regime 

Fully Merged Layer Regime 

(5.1) 


Kn 


*co /R b ■ 


.1 


(3.27) 


and 


e *= .29 


(3.30) 


hence 


V*b « 1 


(5.2) 


So we would expect merged layer flow. Further, we may estimate the shock 
standoff distance by the relation 


« e 

A « 0.5 cm 


(5.3) 


The effect of viscosity would be to substantially increase this value. 
These results correlate with the relatively large freestream Knudsen 
number and Re = 100. 

The effect of ionization or recombination in the flow is to change 
the density of all species. If reactions occur in the shock region 
where the density is changing hydrodynamically , the net result will be 
much more complicated. However, using a simple analysis similar to that 
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used in Chapter III we may justify an assumption of frozen flow through 
the shock. 

Calculation of the ambipolar Schmidt number in Chapter III, 

Sc^ *, 0.6 (3.31) 

indicates that the charged particle density rise at the shock follows the 
neutral hydrodynamic shock. This conclusion may be substantiated further 
by investigating the momentum transfer between ions and neutrals. Colli- 
sion between the heavy particles is effective in momentum transfer since 
the respective masses are nearly equal. Calculating the collision time 
for ion-atom collision: 

T. = — — = 5 x 10 7 sec (5.4) 

ia v ia 


which we compare to the flight time through the shock 


T , = — = 10 sec 

chara v 

s 


(5.5) 


Thus, 


T chara >:> T ia 


(5.6) 


and we conclude that many collisions between heavy particles are suffered 
during passage through the shock. 

It remains to examine the possibility of separation of the electrons 
from the heavy particle shock. We calculate 

\ D « A (5.7) 

from the above results, so this effect is taken as negligible. 
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It is necessary now to investigate the reactions between species. 
Using an analysis similar to that used in Chapter III we may justify an 
assumption of no chemical reactions in the shock. Thus, assuming three- 
body recombination as the dominant reaction, (the two-body recombination 
and ionization coefficients being smaller) we have again from the rate 
equation, 



On 


Az 
e v 


(5.8) 


-12 2 14 

taking a ** 10 cm /sec, N £ = 4 x 10 /cc and Az = 1 cm (~ A) we find 


An e -2 

~ « 4 x 10 (5.9) 

e 


Therefore, the assumption of chemically frozen flow is reasonable. 

To summarize, we expect a viscous shock layer of substantial thick- 
ness. The density rises for the individual species are expected to 
occur at the same position and the flow is frozen with respect to chemical 
composition. From these results we conclude 


N 


N 


e oo 




(5.10) 


Bearing in mind that high electron thermal conductivity may lead to 
elevation of the electron temperature ahead of the shock and consequently, 
an unpredictable temperature increase across the shock, we proceed to 
discuss the optical experiments performed and returned to a more extensive 
analysis of the electron energy equation in the last section of this 
chapter. 
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5.3 Spectroscopic Methods 


The assumption of local thermodynamic equilibrium is frequently used 
in describing the state of a radiating gas. At low densities one may 
usually assume optical transparency (Ref. 47) i.e., at some or all 
frequencies radiation emitted from one part of the gas passes out of the 
gas without being reabsorbed by another part. One then assumes that 
electron collisions with heavy particles and other electrons are suffi- 
ciently dominant over radiative de-excitation that the ionization and 
populations of atomic states are given by the familar Boltzmann- Saha 
equation and that the electrons have a Maxwellian velocity distribution. 

With the recent increase in interest in the spectroscopy of plasmas 
the following more sophisticated theory has developed. The work of Bates, 
Kingston and McWirther (Ref. 39) and others has dealt with the state of 
the plasma when not all the excited states are in a Boltzmann distribution, 
a situation frequently occurring when the electron density is not suffi- 
cient to maintain collisional equilibrium at a given electron temperature. 
The model they assume predicts that when the electrons themselves are in 
a Maxwellian distribution and 3-body recombination dominates, the electron 
collisions tend to maintain population equilibrium of the excited states 
near the ionization limit at the electron temperature. Thus, in the upper 
states , population is determined by collisional effects and the lower 
states are depopulated by radiative de- excitation. The energy level at 
which collisional and radiative transitions become equally probable is 
referred to as the "thermal limit." 

As a criterion for this model, Wilson (Ref. 48) has derived the 
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following relation: if equilibrium is to extend to the u state then 

the electron number density, denoted n fi , must obey the inequality 


n ^ 
e 


r ,*13, w % 7/2 „ m .1/2 
6 x 10 (x - E u ) (kT e ) 

172 


(5.11) 


where x I s the ionization potential, the energy of the u t ^ 1 state and 

kT the electron temperature, all quantities in ev, and n in electrons/cc. 
e e 

It is to be emphasized that this is a necessary but not a sufficient 
condition for upper state equilibrium. Setting = 0 provides an 
approximate criterion for the number densities and temperatures necessary 
for an L.T.E. model. 

Values of n^ are shown for two temperatures in Figure 19. One can 
see that the electron densities needed to maintain equilibrium in all 
excited states are quite high at these temperatures. 

On the basis of this model the most appealing spectroscopic method 
of measuring electron temperature appears to be the Boltzmann plot 
technique applied to the upper states of the atom to obtain the electron 
temperature. This method has been recently reviewed and applied to a 
similar facility in Ref. 49. The authors include much of the existing 
data concerning transition probabilities for argon for convenient 
reference. The principle of the Boltzmann plot is as follows: if a 

transition between state n and m occurs with the emission of radiation, 

1^, at frequency v then 


I = 
V 


g A hv N e 
n nm n o 


E /kT 
n e 


(5.12) 
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where and are respectively the energy and degeneracy of state n; 

tlie trans *- tion probability and N q the neutral species density. 

From Eq. (5.12) if follows that a plot of In ( — ~ — ) vs E^ on semilog 

^n run u 


paper should be a straight line of slope pj--. 

e 

Measurements of electron number density may also be obtained by 
spectroscopy. Again the question of local thermodynamic equilibrium 
becomes important. Calculations based on the Boltzmann- Saha equation 
are ruled out. The intensity of the continuum, however, offers a 
measurement of number density which is at once straightforward and not 
strongly dependent on equilibrium. 

Radiation from the continuum arises from bremstralung (free-free 
transitions) and recombination spectra (free-bound) and consequently 
depend strongly on the velocity distribution of the free electrons. Due 
to the small mass of the electrons and the influence of Coulomb collisions, 
the equilibration time for electrons is very short; therefore, a 
Maxwellian distribution for the electrons in the absence of external 
fields can nearly always be justified (Ref. 50). 

A thorough review of the literature pertaining to continuum radiation 
may be found in Ref. 51. Briefly, the Kramers -Unsold model of the continuum 

predicts that the emission coefficient is given by 

2 

n 


e « 
v 



(5.13) 


where the constant of proportionality is weakly dependent on the frequency 
and the electron temperature. In addition to the assumption of Maxwellian 
electrons the weaker restriction that all of the ions are in the ground 
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state is required. 

Olsen (Ref. 52) has experimentally investigated the argon continuum 
in a free burning arc. For the frequencies between 6000& and 4300& he 
found that the constant of proportionality was relatively frequency 
independent but that the value of the emission coefficient was higher 
than that predicted by the theoretical model. 

Calling the constant of proportionality C and using the relation 
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v 



(5.14) 


we have 
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(5.15) 


This will provide a method of determining the absolute number density 

although the values of C would be expected to be accurate to a factor of 

3. The constant C will be taken from Olsen's data. 

For purposes of studying the shock structure and the region ahead of 

the shock, it is convenient to normalize to the freestream value of n . 

e 

Thus, 



(-V 74 

v Te' 


(5.16) 


This method eliminates the error in determination of the constant of 
proportionality and has the additional advantage of being very insensitive 
to errors in the electron temperature ratio. 
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5.4 The Experimental Arrangement 


For purposes of optical observation the arrangement of equipment was 
as shown in Figure 20. Light from the shock region was folded by two 
mirrors inside the vacuum tank through a viewing port located behind 
the model, then focused at infinity by a 60" focal length lens. This 
lens is mounted on a micrometer drive platform, so that motion of the 
platform corresponds to changing the axial position of the focal point in 
front of the model. The image is then rotated through 90° (for conven- 
ience) and refocused on the entrance slit of the spectroscope, a Hilger- 
Engis grating monochromater (Model 1000, 1 meter focal length) . This 
image was swept across the entrance slits by a mirror which rotates about 
the vertical axis to give a lateral scan. Two different lenses were used 
to refocus the image; one, of 60" focal length gave a 1:1 image, the 
other of 25" focal length gives an image 1:0.24 (the ratio of the entrance 
and exit focal lengths) thus increasing the intensity. The scanning 
speed was varied by moving the final mirror or changing the scanning 
motor. Most scans were about one minute. 

The intensity at the exit slits was detected by an RCA 7850 photo- 
tube mounted in a dry-ice housing to reduce the dark current and powered 
at 2400 V by a Fluke 405B power supply. The output of the phototube 
(voltage) was read directly across the 1 megohm input impedance of a Moseley 
7100A two channel strip chart recorder. The alignment of the optical 
system was accomplished easily by sending a laser beam backwards along the 
optical path starting at the exit slits and ending along the line of 
sight across the body nose. In this manner each of the optics could be 
centered. Determination of the axial position of the body was made at 
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the beginning of a run with the aid of a telescope focused on the 
entrance slits so that the body nose could be centered on the slits. 
Photographs were taken by placing the camera at the position of the first 
lens, so that the same interior optical path was used. 

Because the system is complicated, calibration was performed in 
position. That is, a calibrated tungsten bulb (GE #18A) was placed at 
the nose of the body. Thus, the same optical path is used for both 
calibration and experiments. The lateral scan speed was measured by 
backlighting an object of known size at the flow position and recording 
the output of the spectroscope. The linearity of phototube was measured 
by inserting filters of known optical density into the light from the 
calibrated lamp. The light intensity is decreased in the ratio of the 
power of 10 where D is the optical density. Thus the linearity may 
be measured as shown in Figure 21. The micrometer slit settings were 
checked by plotting the setting as a function of incident to zero 
intensity. The exit slits were found to be accurate within the range 
used. The entrance slits being set much narrower for measurements, were 
generally in error by + lty , scale reading, which was not considered 
excessive. This factor was taken into account during the absolute 
intensity calibration. The calibration given with the standard bulb 
related the brightness temperature to the lamp current. The relation 
between the brightness temperature and true temperature is given by De Vos 
and Rutgers (Ref. 53) who also supply emissivity data for tungsten. We in's 
law is taken to describe the radiation; thus. 
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(5.17) 


ea,T) T 


bulb 


2c 2 h r‘ C 2-, 

"JTT ex P C Tt ] 


8 o o 

where the value of C 2 is 1.438 x 10 A - K . The spectroscope observes lateral 

2 

intensity 1^ (watts/cm - steradian) and the calibration source produces 

2 

intensity (watts/cm - steradian). This information must be converted to 

3 

the radial intensity I (watts/cm - steradian) through inversion and an 

s 

appropriate calibration constant. The inversion is performed in dimensional 
increments of n/R where n is the number of inversion points and R is the 
radius of the jet (3.16 cm). To find the calibration factor a correction 
is needed to account for changes in slit settings between running and calibra- 
tion (Ref. 54). Thus 
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where V = voltage read 

W, W* = exit and entrance slits and 

subscripts c, L denote calibration and lateral signal respectively. The 
inversion scheme then returns 1^. 

For the purposes of Boltzmann plots only the relative intensity is 
necessary. 6500& proved to be a convenient wave length for normalization. 
The response of the bulb was then 
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this quantity was then calculated for several different lamp temperatures 
and applied to the measurements of spectral response described in the 
next section. 


5.5 Spectroscopic Measurements and Results 

Using the apparatus previously described, scans of line intensity 
versus wave length using narrow entrance and exit slits of about 30jj, . 

Some tracings of these are shown in Figure 22. The purpose of these 
scans was to determine line position and width as well as the presence of 
other spectra in the regions of interest. Identification of several 
lines in a given region permitted assignment of wave length units to the 
trace. The appropriate slit setting is then the width of the line (in 
Angstroms) divided by the reciprocal dispersion (8^/mm, first order). 

These results are tabulated in Table II. For the continuum measurements 
the area between 4300X and 4333^ was chosen because of the relative purity, 
thus allowing wide exit slit settings. The cleanest region appeared to be 
centered at 4319$ for a slit setting of 1.5 mm. It was also necessary to 
open the entrance slits to increase intensity. Since this has the effect 
of broadening the line at the exit, these measurements were repeated at 
several entrance settings to insure that the wings of the adjacent lines 
did not overlap the area of interest. Measurements of the continuum were 
conducted with slits of 19C^ width and a 400^ vertical mask. The width of 
the jet was found from the continuum scans to be 6.32 cm. Calibration 
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tests were run at the same time with current and slit settings as follows: 

15 amp, 20jjt entrance, .125 mm exit, V c * 7.5 volts, K = 2.581 x 10"^ 

^ at 2 S 9 Using the measured value of 5.0 x 10 ^ wa ^ ts — at 

cm - ster.volt an ster 


.45" ahead of the body, the emission coefficient 

„ ^-19 watts 

2.7 x 10 — : g j . 

cm - eter - sec* 


e was calculated to be 
V 


From the data of Ref. 49 we calculate 


4319 


7.19 x 10 


■46 watts cm^ 

steradian 


(5.21) 


Using this value and a temperature of 5000°K in Eq. (5.15) gives n^ = 

1.8 x 10 14 /cc. 

For line measurements it was desirable to pick lines which were 
sufficiently isolated to permit opening the slits slightly wider than 
necessary. Thus a small error in wave length setting would not delete 
the wing of the line from the measurements. 

Figure 23 shows a scan of the 4300 line and the underlying continuum. 
The continuum intensity amounts to about 5% of the centerline value of 
the line which was one of the widest measured. Since this error is 
smaller than the reproducibility of a trace, it was decided to neglect 
this correction. In order to retain spatial accuracy and eliminate the 
need for separate calibration of each run, the magnetic field was pulsed 
at approximately one second intervals. Thus a lateral intensity distribu- 
tion with and without magnetic field was obtained, as shown in the 
tracing. Figure 24. This method has the additional advantage of keeping the 
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magnet cooler and the field more uniform. 

In order to recover the radial profiles from the lateral measure- 
ments, the recorder output was first folded and averaged. The method 
used was generally to try to obtain the best agreement about the intensity 
maximum region. Data was rejected if the disagreement in the wings seemed 
excessive. The well-known Abel integral equation was then applied and 
inverted to recover the radial values by a numerical approximation due 
to Pearce (see Ref. 55). Figure 25 shows a sample inversion of some of 
the continuum data; the same technique was applied to the line data. 

For the purpose of making Boltzmann plots, lines were selected which 
has sufficient intensity and isolation and in addition the availability 
of data concerning transition probabilities was taken into account. To 
insure that equilibrium with the freestream could be expected for these 
lines. Figure 19 may be used by placing the upper level energy of the 
line on the graph (as shown for the 4300$ line, the lowest lying line 
selected). The following lines were measured: 4300$, 4333$, 4345$, 4510$, 

5606&, and 5650&. The entrance slits were held to 30^ with both 10C^ and 
40Cp, vertical masks. 

The necessary data for each line used along with the spectral 
sensitivity of the system as measured with the standard lamp is shown in 
Table II. 

The degree of scatter in the final data made construction of a 
straight line difficult, so that a more sophisticated technique was 
required. Accordingly, several sets of data were taken and each set was 
plotted on semi- log paper. Then the data could be superimposed with a 
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relative shift in the direction of the logarithm axis. This corresponds 
to multiplication by an arbitrary constant. Thus, the data were easily 
averaged making construction of a straight line much more accurate. A 
sample Boltzmann plot at 0.15” is shown in Figure 26. The centerline 
temperatures obtained from this method are presented in Figure 27. 
Similarly, ta for the radial variation of freestream temperature was 
shown in re 5. 

Us sse temperature measurements the electron density measure- 

ments may be corrected for the rise in temperature. These results are 
also shown in Figure 26, normalized to the freestream value taken at 
0.35 M . Neglecting the temperature correction Figure 27 shows the variation 
in the shock for several field strengths. While the absolute values are 
not considered highly accurate, as a relative measurement these results 
are considered quite good. The error in this technique is of the same 
magnitude as the scatter in n^ near the centerline (~ 10%) . 

To obtain the absolute values of the electron number density the 
constant of Eq. (5.15) was evaluated. Then with the absolute intensity 
of the continuum known from previous calculations, the electron number 
was obtained. The spectroscopic measurements show expected results in 
that the number density and electron temperature profiles are shifted 
upstream when the field is applied. An unexpected result is that the 
number density increases at the stagnation point. This effect is confined 
to the stagnation region as shown in Figure 28. 
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5.6 Photographic Measurement of Shock Standoff Alteration with Magnetic 

Field 

Photographic techniques were used by Zeimer (Ref. 10) and Seeman 
(Ref. 31) in measuring magnetoaerodynamic shock standoff distances. 

Bailey and Sims (Ref. 56) have presented photographic measurements of 
shock standoff under conditions similar to the flows of this investigation 
at low Reynolds number. Utilizing the self- luminosity of the flow a 
photograph is taken and then scanned with a photodensitometer. They 
assume that the upstream edge of the shock is defined by this point where 
the intensity starts to rise. This is assumed to be the rise in species 
number density. 

A modification of this technique was used to measure the alteration 
to shock standoff. In the region of intensity rise, photographs (such as 
Figure 18) with and without magnetic field show the same initial profile 
(shown in Figure 29). Thus we superpose the two intensity traces with 
and without magnetic field and then move one profile relative to the 
other until they coincide. The distance moved is assumed to be the 
alteration to shock stand-off. This method is illustrated in Figure 
29. 

We are assuming a simple model for the flow in this region, that 
the initial density and temperature rise are being shifted axially but 
do not change shape. This model is reasonable in light of the dis- 
cussion of the preceding section and the spectroscopic data. 


68 



Using this method photographs were taken at various magnetic fields 
at one flow condition. The exposures were made at F/22 and 1/10 sec. 
through a ID neutral density filter. The film was developed for 8 
minutes in Kodak D76 developer. The photodensitometer used was a Joyce- 
Loeble Scanning Microdensitometer. The total system magnification 
accounting for the camera and densitometer was 5.6. The purpose of the 
long time exposure with ID filter was to eliminate some of the intensity 
variations due to plasma noise. 

The sharp intensity rise on the body nose (due to reflection) 
provided a convenient point to measure the alteration. As another measure- 
ment of shock standoff without field, a scratch was placed on the film 
at the nose of the body. The method of Ref. 56 was then followed to 
obtain absolute shock standoff distance. Reasonable agreement was 
obtained. Figure 30 shows the results of this experiment and the altera- 
tion as measured spectroscopically. 

5.7 Observations in the Freestream 

In this section some aspects of the formation of the bow shock wave 
and the partially ionized flow ahead of the shock are considered in more 
detail. In particular, two effects were noticed in the freestream during 
the course of this investigation, a change in the color of the flow (to 
red) when the magnetic field was applied and the "dark space" immediately 
upstream of the shock (see Figure 18) . 

The typical plume of argon plasma produced in this type of facility 
has been discussed by Brewer and McGregor (Ref. 57). The length of these 
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plumes formed as the flow exhausts through the nozzel and into the large 
vacuum chamber is quite long, often several feet, and pinkish-red in 
color. The luminosity of plume is attributed to energy storage in the 
argon metastable excited levels. Since downward transitions are im- 
probable for these levels, atoms excited to these metastable levels tend to 
be convected downstream. However, at any later time electron collisions 
can excite these atoms to higher states from which transitions downward 
to the ground state are possible, thus the energy stored in the metastable 
states is continually lost along the flow. The metastable states of 
argon are at 11.5 ev and 11. 7. The states above these can generate a 
variety of lines through downward transitions, but many of the red lines 
lie within a few ev of the metastable states. The electron temperature 

3 o 

5 x 10 °K) indicates that many electrons available for collisions are 

in this range so that upwards collision to states which have red lines is 
probable, for example, argon 6965$!. In addition, lines such as 6965$ have 
strong transition probabilities so that the argon spectrum in the plume 
tends to be red in color. However, the color change due to magnetic 
field was unexpected and will now be discussed. 

Figure 31 shows several lines of argon in a modified Boltzmann plot 

at 0.65" ahead of the body. Instead of plotting I^/gA as the ordinate, 

Im/I, the ratio of line intensity with and without field is plotted, 

1 T m - T 

giving -r — — — as the slope. Thus, Figure 31 shows an electron tempera- 

K i 1 

m 

ture decrease at this position when the field is applied. The reason for 
the red color change is seen to be a 3- fold increase in the intensity at 
6965^, a relatively strong line of argon in the red. 


70 


N)JW 


To explain the drop in temperature we examine the electron energy 
equation. In one dimensional steady flow following Ref. 41, we have 
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(5.22) 


where 

*= electron thermal conductivity 

* electron viscosity 

• rate of energy transfer by electron ion collisions 

Q « rate of energy transfer by electron atom collisions 
ea 

q *= rate of energy transfer by 3-body recombination 


For the thermal conductivity and viscosity the values for a fully 
ionized gas will be used since the degree of ionization is sufficiently 
high so that Coulomb collisions dominate. Thus for the thermal conduc- 
tivity (Ref. 37) 
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and the viscosity is given by Ref. 38 

4m 

H'e " K e 1st 


(5.24) 


Considering the right side of Eq. (5.22), the rate of electron- ion energy 
transfer has been given by Byron and Petschek (Ref. 59) as 
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(5.25) 


*ei 


". V . 8 ™.,l/2 A ,, ^ »<*V* 
— <TF> ' *> te C ^T? 1 


at m 5000°K, ■ 10^/cc; Q e ^ *= 1,3 x 10^ erg/sec. cm. For purposes 

of calculation, however, we take a simpler approximation (Ref. 37) that 
the rate of temperature change is proportional to the species temperature 
difference and the collision rate 
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dt 
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Thus, 
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We expect that the rate of electron atom energy transfer would be 
a.aaller than the electron ion rate since the respective cross sections 
differ greatly in magnitude while the densities do not. To see this, the 
ratio of energy transfer rates may be expressed approximately as (Ref. 41): 
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where a ea > the atom electron cross section is typically 10 ^ - 10~^cm^ and 

-13 -12 2 

cr the ion electron cross section is typically 10 - 10 cm (Ref, 25). 

ei 

Therefore, we expect Q « Q . and that we may neglect Q in Eq. (5.22). 

ea ex ea 

The three body reaction rate is given by Ref. 44 and Eq. (3.32) as 

dn e 2 

q - € - 5 ! - « % e (3.32) 

where e i8 the energy released during the reaction. The three-body re- 
combination process tends to compete with equilibration since it adds 
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energy to the electron gas. This energy Is released to the third body, 
a free electron during the recombination reaction. Evaluation of the 
reaction rate is made difficult by the possibility that some absorption 
may exist in resonance transitions. This problem has been considered 
by Chen (Ref. 60) who finds that for argon gas at these conditions, 
e " 1 ev, provided the plasma is optically thin, but e can possibly be 
an order of magnitude higher with trapped resonance radiation. We shall 
assume that the plasma is completely optically thin and bear in mind 
that this value of e may be conservative. The value of q is then 

q * 1.7 x 10** erg/sec. cm. 2 


We may further neglect the flow work and viscous terms (containing 
dV 

in the region of interest since V is very nearly constant and in the 
negative z direction. Eq. (5.22) then becomes under the above assump- 
tions: 
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where q find K e will be taken as constant. This equation has an ex- 
ponential homogeneous solution of the form T = Ae"^ 2 where 
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When equilibration is neglected entirely the characteristic length 
for an e-fold decrease in electron temperature 
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10 cm. 


(5.31) 


This is the type of situation considered in (Ref. 33 and 44). Thus, a 
relatively extensive region ahead of the shock would show an elevated 
electron temperature. 
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If equilibration is accounted for, p = 2 cm. from the previously 
estimated values. This value would be somewhat larger if q were taken 
into account. These two values of p = 2 cm with equilibration and 10 cm 
without equilibration represent limiting scale length values for the 
extent of the zone of elevated electron temperature. Since equilibra- 
tion is important as shown previously, the calculations tend to explain 
the fact that such a zone was not observed experimentally (Cf Figure 27) . 

We now consider the interaction of the magnetic field with the con- 
ducting freestream flow. This may be estimated by inserting the Lorentz 
force in a momentum balance. Thus 
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For a field of 200 gauss and a 
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1.5 x 10 mho/meter and Az = 1cm, we 


0.8 (5.33) 


We conclude from this order of magnitude estimate that the flow is slowed 
very slightly by the magnetic field thus forcing the observed temperature 
drop by enhanced electron- ion equilibration. 

Although the flow velocity changes only slightly the equilibration 
rate is sufficient to cause a change in the relative line populations. 

This idea is further substantiated by Figure 32 which shows that the 
6965$ line of argon (which shows the strongest intensity change with 
magnetic field) is very strongly altered even at the outer fringes of the 
jet. This would be expected if one considers that at the outer regions the 
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magnetic field lines are more nearly normal to the flow, thus the flow 
would be even more slowed here than at the centerline. 

The cause of the "dark space" in the luminosity just ahead of the 
shock observed here and in other investigations then remains to be 
explained. Grewal and Talbot (Ref. 33) attribute this to the zone of 
elevated temperature causing a decrease in recombination, hence in line 
radiation from the plasma. 

If we assume as they do that the 3-body recombination coefficient 
is given by 

a = 5.6 x lo" 27 (Te) _9/2 n e (5.34) 

then the intensity ratio of freestream to dark space is given by 

T 

i I -i CC 9 - eoo /c 

log 10 I " log 10 a “ 2 log 10 Te (5.35) 
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The value of log I/I may be found from Figure 29, since 

00 

log 10 f " D - D w (5.36) 

00 

the film density change. Thus we find 
T 

T = 1.1 (5.37) 

00 

so that an extremely small rise (10%) in electron temperature while the 
number density remains constant can account for the "dark space" in the 
manner offered by Ref. 33. This type of accuracy is beyond the accuracy of 
Boltzmann plot technique we have used, so that the explanation can neither 
be confirmed or discounted. However, the fact that such a small difference 
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which are then solved for k^. This method is compatible with the concept 
of standard deviation of the data which is then defined (Ref. 61) 


6 



( 6 . 2 ) 


where n is the number of discrete data points. 

By inspection of the data and because no directly applicable expres- 
sions exist, the data for the force on the magnet and the pressure altera- 
tion are fitted with linear relationships. Since considered as a percentage 
alteration these quantities must pass through the origin, the following 
function is taken. 

f = aB 2 (6.3) 


Applying the above process results in 

Lorentz drag force on magnet 
Total no field drag on body 


x 100% = 1 . 88B 2 % (6.4) 


Alteration to Pressure force inn7 _ , Q 2 0/ 
Total no field drag on body X ° ~ " * 


(6.5) 


Inspection of the data (Figure 13) also shows that it is not reason- 
able to apply a linear function to the total drag alteration. Since the 

2 

Hall effect (which depends on B rather than B ) is expected to be import- 
ant a series of the form 

f <= ajB + a 2 B 2 + a 3 B 3 + (6.6) 


could be considered. A two term expansion of this form however, proved 
to give a bad fit, although the general features were reproduced. While 
higher order expressions chould be taken, physical interpretations of 
empirical expressions derived in this manner are difficult to make. A 
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simple function of the form 


f» — T (6.7) 

1 + a x 

may be fitted to the data. The motivation for picking this particular 
form will be discussed later in this section. The following procedure was 
used to fit the function. The parameter may be determined graphically 
by trial and error since the function is less sensitive to than to a; 
f is then linear in X where 


X = 


1 + a, B" 


( 6 . 8 ) 


The constant a may readily be evaluated as before. The function f deter- 
mined is 


2 

alteration to total drag infvy __ 1.55B ^ 0/ 

total no field drag on body x ° 1 + .01B ° 


(6.9) 


The data for the alteration to the shock standoff as 

30 was also fitted to a linear relationship. Thus, 

alteration to shock standoff 

■ — = 0.96 x 10 

body radius 


shown in Figure 
B 2 % (6.10) 


The relationships for the drag and shock standoff alteration may be 
expressed in terms of the interaction parameter given by Eq. (3.30). The 
resulting expressions along with the standard deviation 6 about the curve 
fit are summarized in Table 3. 

These expressions are shown plotted with the data in Figure 33. In 
addition to these functions, we may deduce a function for the viscous drag 
alteration by subtracting the measured components from the total drag 
alterations to the measured components. 
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Thus 



This would be reasonable for the range of B =10 to 20, where the func- 
tion is negative. 

The function (Eq. 6.9) selected to represent the total drag was some- 
what unusual and should be discussed further with regard to its behavior. 
Judging from the calculated standard deviation this function provides a 
good fit, which is sufficient to justify using it to represent the data. 
However, this function is also attractive since some theoretical arguments 
have been put forth to indicate that beyond a certain value (shock lift- 
off) the value of the drag can no longer increase due to the fact that 
the magnetic force is balancing the force of the flow through the shock 
(Refs. 3 and 25) „ While it is generally recognized that this idea is 
largely hypothetical based on elementary considerations and does not account 
for such effects as change in shock shape, there is some validity to the 
argument especially in the laboratory where the flow is not unbounded and 
there is a maximum value for the flow momentum. We would expect then 
that the increase in drag would taper off at high interaction and possibly 

approach a constant value. The function chosen for the total drag data 

2 

has this behavior since it tends to a constant value as B becomes large. 
This value of the function represents the point at which shock layer lift- 
off could conceivably occur and no further increase in drag may be pos- 
sible. The ultimate value of the function is 155%. We may estimate the 
interaction necessary to produce this type of drag alteration and the 
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approach to shock lift-off. For example, at S = 200, the drag alteration 
would be one-half its final value and at S = 400, two- thirds. 

6.2 Comparison of Magnetoaerodynamic Drag Data to Theory 

The results in Ref. 25 are presented as magnetoaerodynamic drag 
coefficients which are functions of the interaction parameter S. These 
values may easily be converted to a percent alteration using the value 
of the coefficient at S - 0. The empirical curve (Eqs. 6.9 and 3.30) for 
total drag is shown in Figure 34 with the following results from the 
theory of Ref. 25: inviscid, no Hall effect; inviscid, with Hall effect; 

viscous, no Hall effect; and a composite theory accounting for both the 
Hall and viscous effects at high interaction. In order to correctly refer 
the data to these theories it was necessary to account for the difference 
between e = .29 for experiment and e = .20 for the theories. This was 
done by correlating the data and theory by the nondimens ional parameter 
Se as explained in Chapter I. Thus, the different theories could be 
readily converted to e = .29. In the Hall effect theories the Hall 
coefficient varies with B in the same manner as experiment. 

The composite theory of Ref. 25 is based on the theoretical result 
that the viscous solution merges with the inviscid solution for moderate 
interaction. It is supposed that this also occurs when the Hall effect is 
present, so that for moderate interaction the inviscid. Hall effect drag 
coefficient is taken to be the solution. To construct the percent 
alteration then the following expression is taken: 
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x 100% (6 . 15) 
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Since C (0,0,R ) is larger than C (0,0,R = oo) the theory will show a nega- 

1 / e i/ e 

tive alteration if S is taken too small. Thus, the composite theory does 
not pass through the origin in Figure 34. 

Inspection of the viscous and composite theories in Figure 34 shows 
that the onset of substantial increases in total drag is considerably 
delayed by viscosity. This theoretical result has been attributed to a 
near cancelation of effects (Ref. 25). The pressure and viscous drag 
coefficients both increase as the Reynolds number decreases in purely 
aerodynamic flow. At a particular Reynolds number the effect of magnetic 
interaction is to cause an increase in the Lorentz force which is nearly 
cancelled by a decrease in pressure and viscous drag. Presumably 

the viscous coefficient goes to a constant value at some value of S and 
so can cause no further effect while the force on the magnet keeps 
increasing. This explains qualitatively the delayed feature of the viscous 
theory curve. 

While experimental conditions indicate that the flow is very viscous, 
the experimental drag data does not reproduce the qualitative feature of 
the viscous theory curve of Porter (Ref. 25). The best qualitative agree- 
ment is found with the inviscid theories, and at high interaction parameter 
with the composite theory. For the present we will not be concerned with 
quantitative agreement and try to understand more completely the qualita- 
tive features of the data. 

In the experiment it was found that the viscous drag on the body was 
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changed bery little when the magnetic field was applied. Eq. (6.11) 

2 

would show that at 25 Kg the alteration to the total drag as a result of 
alteration to the viscous drag is about 4 percent. This can be explained 
by the fact that the axial component of the viscous force near the stagna- 
tion region is small as compared to the total axial viscous force for an 
hemispherical body with an afterbody. This is because the axial component 
of the shear stress at the stagnation region is small and the stagnation 
area is likewise small. Since the electromagnetic effect is concentrated 
only near the stagnation region, it is perhaps not surprising to see that 
the change in viscous drag due to magnetic field is small as compared to 
total viscous drag. 

According to the experimental evaluation of the interaction para- 
meter, a 4 percent change in total drag due to viscous effect would cor- 
respond to an interaction parameter of 50. Viscous theory however, pre- 
dicts that the viscous drag would have approached a constant value at 
this interaction parameter since this is the range where the inviscid 
and viscous drag coefficient are nearly identical and the viscous theory 
predicts a substantial increase in drag with interaction. This is also 
the range where the composite theory may begin to be applicable. 

Part of the reason for the discrepancy between the viscous theory 
and experiment is due to the small angle assumption of the theory. In 
the region near the stagnation point, the assumption of sin 9 0 is 

justifiable. The extension of the calculation to 45° using the above 
assumption would tend to overestimate the change in viscous drag due to 
magnetic field. Another factor that influences the theoretical result is 
that Porter only calculated viscous drag up to 45°, and neglected the 
viscous drag at large angles and of the afterbody. This would tend to 
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reduce the total drag coefficient thus dramatize the decrease in viscous 
drag coefficient. 

Still another factor may contribute to the discrepancy and that is 
the assumption of a thin shock in the theory. At very low Reynolds 
numbers encountered in the experiment the shock is actually quite thick. 
Thus, the actual flow profile around the body may be quite different from 
a theoretical profile calculated assuming a thin shock. 

These reasons indicate why the experiment does not have the quali- 
tative feature of the viscous theory curve in Figure 34. The viscous theory 
over-estimates the amount that the viscous component of drag is altered 
when the magnetic field is applied. While the velocity profile around 
the body is certainly different than the inviscid theory would predict, 
the principle alteration of the drag on the body is the Lorentz force 
and the alteration in pressure force. The viscous force stays relatively 
constant because it is important mostly outside the region of interaction 
around the nose. The reason that the inviscid theory and the data do 
not agree quantitatively is that the theory which gives the alteration to 
drag as a percentage is based on the inviscid drag coefficient at no 
magnetic field whereas the experiment is based on a viscous drag coeffi- 
cient. Other reasons are that the velocity profiles differ and the experi- 
mental determination of S 'may be off. Suitable modifications to the 
theory to account for quantitative discrepancies will be discussed later. 

We now discuss comparison of the data with theories which include 
the Hall effect. We have noted previously that the total drag data curve 
agrees qualitatively with the inviscid Hall effect theory and also with 
the composite theory at. high interaction. However, to fully understand 
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whether or not this behavior is completely in agreement with the theory 
it is necessary to examine individually the Lorentz and pressure forces. 

We consider first the Lorentz force. In Porter's theory (Ref. 25) 
this force is found to be independent of viscosity. Apparently this 
result is again a cancelation of effects. The velocity for the more 
viscous cases is compensated by the increase in the interaction volume 
due to the increase in shock standoff with viscosity. Because of this 
theoretical result it is easy to formulate the percent alteration to 
the total drag for the Lorentz force component including both the Hall 
and the viscous effects. This is done in a manner analogous to the com- 
posite theory for the total drag alteration. Thus, 


% 


alteration, Lorentz component = 


C DL^ S,C H» R e ~ 00 ) 
C D (S=°, C H =0,R e ) 


x 100% 


(6.16) 


Figure 35 shows the empirical relation for the alteration of the total 
drag force due to the force on the magnet plotted with the theory for 
C u = 0 and C„ = .56 JS y the value of C„ determined experimentally. The 
reaction force on the magnet is found to increase more rapidly with S 
than would be predicted by the Hall effect theory. 

The pressure component of total drag may be examined in a similar 
manner. Figure 36 shows the experimental curve for the alteration in 
total drag due to change in pressure drag and three theoretical curves, 
a viscous and inviscid theory (no Hall effect) and an inviscid Hall 
effect theory. The manner in which the theoretical curves are generated 
was explained in Chapter II. 

Figure 36 shows a surprising qualitative difference. The inviscid _ 
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Hall effect theory shows that the pressure force on the body increases 
with interaction whereas the non-Hall effect theories show a decrease. 

The experimental data shows that the pressure force decreases. 

In addition to the pressure force on the body, we may also examine 
the experimental pressure distributions. (Cf. Figure 17). These show 
that the pressure is locally relieved at angular positions where the 
projected area is large (^20° to 40°). The stagnation pressure rise 
contributes only a slight increase in pressure force since the projected 
area over which this pressure acts is small. 

We cannot explain why the inviscid Hall effect theory of Ref. 25 
shows that the pressure force on the body is increased. However, it 
may be noted that this does not seem to be compatible with the notion 
that the shock layer may eventually lift completely off the body and 
totally relieve the pressure. 

In addition to the disagreement between theory and experiment, with 
regard to the pressure force, the observed pressure rise at the stagna- 
tion point would not be expected from theoretical results since it is 
assumed that there is no interaction along the central streamline. A 
plausible explanation of this phenomenon is that electrons tend to follow 
the magnetic field lines. Thus, for convergent field lines, as in the 
present case of a dipole magnet at the center of a hemisphere, the elec- 
trons are being channeled toweard the stagnation region. At the same 
time the electrons also follow the fluid motion along the streamline 
since we have charge neutrality and no ion slip. In a weak magnetic 

field with C„ < 1, the collisions between electrons and ions and neutrals 
n. 

are strong enough so that the electrons move along the streamlines. In a 
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strong magnetic field with > 1, as in the present case, the drifting 
of the electrons along the field lines becomes importnat, Because of 
charge neutrality and no ion slip, the ions and neutrals also tend to be 
pulled toward the stagnation region* The net result is an increase in 
density in the stagnation region resulting in an increase in temperature 
and pressure in the stagnation region* 

From this and the previous discussions we conclude that while some 
agreement seems evident between the various theories and experimental 
data for total drag, closer inspection of the various components of the 
drag has shown that this agreement is not good qualitatively. 

The theoretical and experimental results for shock standoff altera- 
tion are compared in Figure 37* Here the inviscid and viscous theories 
lie close together and are qualitatively similar to the data. The ex- 
perimentally measured alteration to shock standoff seems to be less than 
that of either the viscous or inviscid theories. 

Another reason for this difference is undoubtedly the influence of 
the Hall effect which is accounted for in both the inviscid and composite 
theories. However, we did not observe the inflection in the data which 
both these theories predict. Again no reason has been given for these 
inflections. Possibly the scatter in the data was too large and magnetic 
field too low to permit this observation. Because the shock structure in 
experiment and theory differs substantially it is difficult to make a 
meaningful quantitative comparison. It would seem that the non-Hall 
effect theories predict too high an alteration while the Hall effect 
theories predict too low an alteration. No arguments against the viscous 
theory similar to those made in the discussion of the drag can be made 
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here. However, the Hall effect theory which is already under question 
concerning the drag should also be questioned here because of the unex- 
plained inflection in the alteration. 
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CHAPTER VII 


SUMMARY AND CONCLUSIONS 

In the course of this investigation, we have observed the phenomena 
of magnet oaerodynamic flow for a simple blunt body problem. Observation 
of the expected increases in magnetoaerodynamic drag and shock standoff 
has been combined with supplementary investigation of the non- equilibrium 
properties of the ionized gas flow. 

Spectroscopic measurements in the shock layer and freestream have 
been made. The fact that the flow around the body is very viscous was 
confirmed by experimental measurements of the shock profile. This pro- 
file was found to be similar to theoretical models for merged layer flow 
(Cf. Ref. 62). In the freestream and in the shock the flow was non- 
equilibrium. While the drag is not affected by the non-equilibrium free- 
stream, the small amount of interaction is seen to account for the ob- 
served color change when the magnetic field is applied. It should be 
emphasized that this color change was observed in argon and depends on 
the metastable states of argon. Thus, for a different gas other color 
changes or no color change at all may be observed. Behind the shock the 
non-equilibrium effects tend to keep the electron temperature constant. 
Thus, the assumption of constant conductivity seems quite reasonable for 
this experiment. 

The experimental drag measurements show an substantial increase in 
total drag at a magnetic field of several thousand gauss. Further meas- 
urements show that the increase is due to the increase in Lorentz drag, 
with the pressure drag actually decreasing and the viscous drag decreas- 
ing slightly. The shock standoff shows a substantial increase with 
magnetic field. 
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The pressure and spectroscopic measurement furthermore show that 
pressure, density and temperature of the gas increase with magnetic 
field near the stagnation region. This is analogous to "mirror effect" 
and indicates that Hall effect is important in this experiment. 

An attempt is made to compare the experimental result with different 
calculations especially that of Porter's (Ref. 25). Although both the 
calculation and the data agree in overall trend of increase in drag and 
shock standoff with magnetic field, they show some disagreement in de- 
tails. They are: 

1. Experiment shows an increase in density and pressure near the 
stagnation region. The constant density model by its assump- 
tion would exclude such changes. 

2. Experiment shows a decrease in pressure drag with magnetic 
field and theoretical calculation with Hall effect actually 
shows an increase with magnetic field. It is not known at 
present why the theoretical calculation shows an increase of 
pressure force with Hall effect and yet a decrease when Hall 
effect is neglected. 

3. In the shock standoff measurement, again the data do not agree 
qualitatively with the Hall effect theory. No physical ex- 
planation has been given for the existence of the inflection 
point on the theoretical curve. 

In conclusion the present experimental investigation shows that 
substantial increase in drag with high magnetic field is possible. The 
experiment further shows that prediction based on present theories should 
be treated with care. For high magnetic field with Hall effect , constant 
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property model breaks down. As yet there is no completely satisfactory 
Hall effect theory. For even higher magnetic field to achieve lift-off 
of the shock layer, depending on the degree of ionization, ion slip can 
play an important role. Thus any conplete theory should take this into 
account at high magnetic field as in Hall effect in moderate magnetic 
field. 
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Pig. 2 Schematic of facility 



Fig. 3A Photograph of facility showing operating console 
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Fig. 4 Photograph of pressure probe and various heads 


103 


IMPACT PRESSURE (mm Hg) 


o N e (.35"AXIAL) 

o FREE STREAM IMPACT 
PRESSURE 



Fig. 5 Freestream property profiles 
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Fig. 6 Interior detail body and magnet 





TOTAL DRAG MODEL 
AND MAGNET 


I 




FLAT END MOOEL 
AND MAGNET 



PRESSURE MODEL 


I 



LORENTZ DRAG 
MODEL 


Fig. 7 Photograph of bodies and magnets used in magnet oaerodynamic 

drag measurements 
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Fig. 8 Photograph of drag measuring apparatus 
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Fig. 10 Sample magnetic field strengths 
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Fig. 11 Field drop off with axial distance 
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Fig. 13 Alteration to total drag and components with magnetic field 
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Fig. 14 Comparison of total drag alteration for two body 
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Fig. 15 Pressure distribution on body for two tap sizes 




Fig. 16 Pressure change at various angles as a function of magneti< 
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Fig. 17 Pressure distribution on body for several field strengths 
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Fig. 18 Photograph of body in flow with and without magnetic field 
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Pig. 21 Response of phototube 


120 










INTENSITY (ARBITRARY UNITS) 



NONDIMENSIONAL JET RADIUS 

Fig. 25 Typical inversion of lateral scan 
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Fig. 26 Sample Boltzmann plot 
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Fig. 28 Radial profiles of electron temperature and number density 
in the stagnation region 
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Fig. 29 Showing the photographic technique used to measure 
shock standoff alteration. 
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Fig. 31 Relative Boltzmann plot at 0.65" axial position 
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Fig. 34 Comparison of total drag data and theory 
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Fig. 35 Comparison of Lorentz drag data and theory 
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PERCENT ALTERATION TO TOTAL DRAG 



Fig. 36 Comparison of theoretical and experimental pressure force 

measurement 
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Fig. 37 Comparison of shock standoff data and theory 



APPENDIX A 


Facility 

A detailed schematic diagram of the facility is shown in Figure Al# 
The specifications for the equipment are as follows; 

Vacuum System 

Main pump: Roots Connersville 

Rotary gas pump (impeller) 

Type RGS 
Size 10 x 30 
5,000 cfm 

Backup pump: Kinney Corporation 

Rotary pump (positive displacement) 

Model KD485 
460 cfm 

Ambient pressure gauge: 

"Alphatron" 

National Research Corporation 
Type 530 
(Matching head) 

Vacuum line: 8" diameter steel pipe throughout 

Neophrene gasket seals 

Main valve manufactured by Hills McCanna Corporation 
(Bonnet Type) 
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Arc Heater and Power Supply 


Plasma Torch: Thermal Dynamics Corporation 

Model F80 

(Maximum power input 80 KW) 

Standard electrode and anode 91B, 1/2" diameter 
Power Supply: Miller 

(welders) 

Model S.R.H. .444 
(4 units) 

75 V.D.C., 400 amps, each 
Gas Flowmeter: Fisher and Porter 

Flowrator 
Model 1 0A 3500 

Some typical operation characteristics of the facility are shown in 
Figure A2 as a function of argon flow rate for a current setting of 960 
amps. While a substantial change in ambient pressure could be caused by 
shutting off the blower pump the jet obtained was too small to be useful 
(~ 1” diameter) . 

Auxilliary equipment available for use with the facility included 
the following: 

1. Davenport Power Supply, Model VD2-30K-1E, 30 amp/200 VDC (used 
for magnet power and battery charging) . 

2. Volta Batteries - 60 amp hr. (storage type) -- 12v (32 used as 
main magnet power supply.) 

3. Harrison Labs Model 6264A, 0-30 amp DC (constant current- constant 
voltage used for optical calibration). 
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4. Fluke High Voltage Power Supply Model 405B, 0-2700 v/30 ma 

(Photomultiplier power supply). 

5* Hewlett Packard/Dymec Scanning Digital Voltmeter: 

Model 2401B Integrating Voltmeter 
Model 2901A Input Scanner 
Model 562A Digital Recorder 

6. Moseley Two Channel Strip Chart Recorder, Type 7100A. 

7. Bell Gaussmeter Model 120 (0-30 KG) 

The drag measurement system described in Chapter IV is shown 
schematically in Figure A3. The equipment specifications were as follows:*" 

1. Linear Differential Transformer: 

Sanborn Model FTA 100 
Capacity 0-100 grams 
Sensitivity 8 mv/volt applied 
Linearity 0.2% fullscale 

2. Power Supply and Demodulator: 

Sanborn Transducer Converter Model 592 
(operating at 2.4 KC/5 volt DC input) 

3. Recorder: 

Moseley 7100A two channel strip chart recorder. 


1. The optional filter circuit used to eliminate mechanical noise is 
described in "Andreyev's Circuit," Electronic Engineering , Vol. 37, 
p. 444, February, 1965. 
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Fig. A1 Schematic of facility 
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Fig. A3 Schematic of drag measuring device 
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APPENDIX B 


Diagnostic Measurement 

The method of determining impact and static pressure was discussed in 
Chapter III and the probes used for these measurements are shown in 
Figure 4. The interior detail of the sting support is shown in Figure 
Bl. 

The pressure sensor was a daiphram type, the displacement of the 
diaphram being read by a linear differential transformer mechanism. The 
specifications are as follows: 

Sanborn Model 270 Differential Pressure Transducer: 

Linearity 1% full scale 
Sensitivity 28 MV/volt applied 
Minimum Differential Pressure 1.5 ^ 

Sanborn Transducer Converter Model 592: 

(Operating at 2.4 KC/5 volt DC input) 

Some typical radial surveys made with the impact pressure probe under 
three different conditions are shown in Figure B2 . 

The mass flux probe tip is shown in Figure B3 along with the schematic 
of the experimental arrangement. To measure the pressure in the exterior 
dump tank the following "Alphatron" gauge was used, 

"Alpha tr on" vacuum gauge: 

National Research Corporation 
Model 520 


(matching head) 



The spectroscope used was: 

Hilger-Engis Model 1000 
1 meter focal length 

(Grating type monochrometer 5000j? blaze-1200 lines/mm (8$ first 
order)) 
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Fig. B1 Impact pressure probe 
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Fig. B2 Typical Impact pressure surveys 
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Fig. B3 Mass flux experiment -- detail. 
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